UNCLASSIFIED 


f 


Reproduced 
luf.  the 


ARMED  SERVICES  TECHNICAL  INFORMATION  AGENCY 
ARLINGTON  HALL  STATION 
ARUNGTON  12,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever]  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


n 


[T4^  »“ 

'  i 

j 

PRINCETON  UNIVERSITY 

DEPARTMENT  OF  AERONAUTICAL  ENGINEERING 


U.S.  Army  Transportation  Research  Command 
Fort  Eustis,  Virginia 

Project  Number:  9-38-01-000,  TK902 

Contract  Number:  DA  44- 177-TC-524 


STABILITY  AUGMENTATION 
OF  GROUND  EFFECT  MACHINES 

by 

Theodor  A.  Dukes 
Charles  R.  Hargraves 

Department  of  Aeronautical  Engineering 
Princeton  University 
Instrumentation  and  Control  Laboratory 

Report  No-  601  April  1962 


FOREWORD 


The  research  in  this  report  was  conducted  in  the  Instrumentation 
and  Control  Laboratory  of  the  Department  of  Aeronautical  Engineering 
at  Princeton  University. 

The  work  was  done  under  the  sponsorship  of  the  U.S.  Army  Trans¬ 
portation  Research  Command  as  Phase  8  of  ALART.  This  support  is 
gratefully  acknowledged. 


TABUS  OF  CONTENTS 


page 


LIST  OF  FIGURES . 1 

LIST  OF  SYMBOLS . 3 

SUMMARY . 7 

1.  INTRODUCTION  .  8 

1.1.  Background  . . 8 

1.2.  Equations  of  motion . 11 

2.  HOVERING  ANALYSIS . l4 

2.1.  Hovering  equations . l4 

2 .2 .  Application  of  momentum  balance  to  GEMs . «  .  15 

2.2.1.  Assumptions . 15 

2.2.2.  Balanced  regime  . . l6 

2.2.3.  Underfed  regime . 18 

2.2.4.  Overfed  regime  .  18 

2.3.  Lift  and  augmentation . 19 

2.4.  Heaving  dynamics  .  21 

2.4.1.  Restoring  force  in  heave  .  21 

2.4.2.  Damping  in  heave . 24 

2. 4. 2.1.  Downward  velocity  . . 24 

2. 4.2 .2.  Upward  velocity  . . 26 

2.4.3.  Frequency  and  damping  ratio  of  the  heaving  motion  .  .  29 

2.5.  Pitching  motion . 30 

2.5.1.  Moments . 30 

2.5.2.  Three  dimensional  correction  .  42 

2.6.  Determination  of  base  pressure  and  jet  momentum . 45 


page 


3-  DISCUSSION  OF  THE  THEORETICAL  AND  EXPERIMENTAL  RESULTS . 49 

3.1.  Heaving  motion . 49 

3.2.  Pitching  motion . 51 

4.  FORWARD  FLIGHT  CONSIDERATIONS . 57 

5.  PRELIMINARY  CONTROL  SYNTHESIS . 66 

5.1.  Stability  augmentation . 66 

5.1.1.  Description  of  the  uncontrolled  system  .  66 

5.1.2.  Design  principles  for  the  controlled  system  .  69 

5.2.  Preliminary  design  of  a  feedback  system  for  attitude 

stability  augmentation  .  74 

6.  CONCLUSIONS . . . 86 

REFERENCES . 89 

APPENDIX  A . 92 

APPENDIX  B . 96 


FIGURES 


LIST  OF  FIGURES 


1.  Comparison  of  Jet  Theories 

2 .  Lift  versus  Height 

3.  Equilibrium  Height  versus  Attitude  Angle 

4.  Heave  Frequency,  0°  Jet  Inclination 
5-  Heave  Frequency,  45°  Jet  Inclination 

6.  Experimental  Heave  Frequency 


7 .  Heave  Damping 

8.  Experimental  Heave  Damping,  =  3.86 

9.  Experimental  Heave  Damping,  =2.02 

10.  Experimental  Heave  Damping,  R_  =4.78 

11.  Experimental  Heave  Damping,  4350  RPM 


12 .  Heave  Mode  Root  Locus 


13.  Effect  of  Cross  Flow  Dissipation  Factor  on  >  Hu  =  0.7 

14 .  Effect  of  Cross  Flow  Dissipation  Factor  on  ,  Ql.  =0.8 


15.  Effect  of  Nozzle  to  Base  Area  Rotor  on 


dCy 

ae 


l6.  Effect  of  Jet  Inclination  Angle  on  ~ 


5>Cm 

Z>e 


1 


17.  Effect  of  Cross  Flow  Dissipation  Factor  on 

18.  Effect  of  Lift  Efficiency  Factor  on  -L  ,  P,o  =  0° 

*  ©  ©  ' 

19.  Effect  of  Lift  Efficiency  Factor  on  -£■  =  ^5° 

20.  Experimental  Values  for  -L  ***&& 

21.  Experimental  Moment  versus  Attitude  Curves 

22.  Experimental  Pitch  Damping 

23.  Change  of  the  Pitching  Mode  with  Height 

2k.  Effect  of  RPM  on  the  Pitching  Mode  (Experimental) 

25.  Effect  of  Damping  on  a  Statically  Unstable  Second  Order  System 

26.  Feedback  Configurations 

27.  Close-up  View  of  the  Model 

28.  Model  and  Test  Setup 


2 


LIST  OF  SYMBOLS 


A 

Ac 

af 

Ap 

Ar 


Oo 

B 

Cj>p 

CL 

CM 

cto 

Cv, 

D 


Vi 


gain,  augmentation  factor 

control  gain 

attitude  feedback  gain 

pilot  gain 

rate  feedback  gain 

nozzle  area  to  base  area  ratio 

damping  force  ( lbs ./ft ./sec . ) 

parasite  drag  coefficient 

lift  coefficient 

moment  coefficient 

momentum  coefficient 

factors  defined  by  Equation  22 
equivalent  diameter  (ft.) 
acceleration  of  gravity  (ft ./sec.2) 
general  transfer  function 
height  above  the  ground  (ft.) 


3 


\o 

I 

J 

I 

K 

L 


m 

P 

P 

Pu 

r 


normalized  height  =  — 


moment  of  inertia  (slug-ft.2) 
jet  momentum  (lbs. /ft.) 

T 

normalized  jet  momentum 


AfVt 

aerodynamic  spring  constant  (lbs. /ft.) 
lift  force  (lbs.) 


length  of  the  periphery  (ft.) 

M  aerodynamic  moment  (ft. -lbs.);  mass  of  GEM  (slugs) 


m  mass  flow  ( slug/ft . -sec . ) 


normalized  mass  flow  = 


pressure  (lbs. /ft.2) 


normalized  pressure  = 


specific  load  =  -=r  (lbs. /ft.2) 


normalized  specific  load 
radius  or  distance  (ft.) 


radius  of  gyration  of  saucer  (ft.) 


V"1  normalized  radius  = 


3  base  area  (ft.2) 


■t 


k 


S  complex  variable 

SC  control  lag  root 

S,  ,  S2.  pitching  mode  roots 

T  propulsive  thrust  (lbs.);  time  constant  (sec.) 

^  jet  thickness  (ft.) 

V  velocity  (translational)  (ft ./sec.) 

(VT  jet  velocity  (ft. /sec.) 

w  weight  of  GEM  (lbs.) 

a/  _  \  +  Awn  £ 


o< 

p 

S 

S' 

It 

n 

© 


cross  flow  dissipation  factor 

jet  inclination  angle  with  respect  to  the  vertical  (deg.) 

jet  inclination  angle  with  respect  to  the  machine  (deg.) 

control  input 

damping  ratio 

jet  efficiency  factor 

base  (lift)  efficiency  factor 

attitude  angle 


5 


_  - 

— 

| 

1 

i 
| : 

1  j 

. 

p 

density  of  air  (slug/ft. 3) 

1 

.5 

1 

CO 

frequency  ( 1/ sec . ) 

1 

SUBSCRIPTS 

A 

atmospheric 

B 

base 

b 

balanced  regime 

J 

e 

equilibrium 

M 

high  side  or  heaving  mode 

« 

T 

jet 

L 

low  side 

o 

overfed  regime,  free  stream 

P 

pitching  mode 

T 

total;  throttle 

u. 

underfed  regime 

oO 

out  of  ground  effect 

© 

angle 

1 

6 

% 

SUMMARY 


This  is  a  study  of  feedback  control  for  the  stabilization  of  a  peri¬ 
pheral  jet  ground  effect  machine  (GEM).  The  study  is  limited  to  over  land 
operation,  to  hovering  and  low  forward  velocities,  where  normal  aerodynamic 
control  surfaces  are  ineffective. 

Expressions  for  the  frequency  and  damping  of  the  heave  motion  and  for 
the  attitude  moment  derivative  are  derived,  using  the  principle  of  momentum 
balance.  The  influence  of  physical  parameters  and  scaling  on  these  results 
is  discussed.  The  differential  equations  of  forward  flight  are  developed 
and  expressed. 

The  results  of  a  series  of  experiments  on  an  eight  foot  diameter  GEM 
model  in  a  hover  condition  are  in  satisfactory  agreement  with  theoretical 
predictions.  The  experiments  also  show  considerable  damping  of  the  attitude 
motion. 

An  attitude  and  rate  feedback  control  system  was  devised,  considering 
the  adaptability  of  the  human  pilot,  and  the  disturbing  moments.  Pre¬ 
liminary  synthesis  brings  out  the  significance  of  the  moment  control  lag 
and  the  variation  of  the  moment  control  effectiveness.  It  is  suggested 
that  open  loop  gain  adjustments  can  be  expected  to  provide  satisfactory 
compensation  for  parameter  changes,  both  in  hovering  and  in  forward  flight. 
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INTRODUCTION 


1  . 

1.1.  Background 

Most  low  speed  flying  vehicles  are  either  unstable  or  exhibit  only 
marginal  stability.  Two  outstanding  examples  of  vehicles  with  stability 
problems  are  ground  effect  machines  (GEMs)  and  vertical  take  off  and  land¬ 
ing  aircraft  (VTOLs).  Although  it  is  a  well-known  fact  that  aircraft  need 
not  be  actually  stable  (i.e.  return  to  equilibrium  after  any  disturbance), 
for  satisfactory  flying  qualities  it  is  necessary  that  any  divergence  in 
the  motion  following  a  disturbance  be  sufficiently  slow  so  that  it  can  be 
comfortably  handled  by  the  pilot.  Many  GEM  and  VTOL  configurations  are  suf¬ 
ficiently  unstable,  over  at  least  a  portion  of  the  flight  regime,  as  to  be 
difficult  if  not  impossible  to  fly.  Aerodynamic  stabilization  in  many 
cases  proves  to  be  inefficient.  Another  means  for  increasing  the  stability 
is  the  application  of  feedback  control.  Some  difficulty  may,  however,  be 
encountered  in  the  application  of  conventional  techniques  to  automatic  stabi¬ 
lization  because  of  the  rapid  variation  of  the  stability  derivatives  during 
changes  in  flight  condition. 

The  report  which  follows  is  confined  to  the  problems  of  artificial  sta¬ 
bility  augmentation  of  peripheral  jet  ground  effect  machines.  An  investigation 
of  the  literature  shows  that  very  little  is  known  about  the  dynamic  character¬ 
istics  of  this  type  of  vehicle. 

Attitude  stability  can  be  obtained  in  the  GEM  by  aerodynamic  means .  If 
the  GEM  is  operated  at  a  height  above  the  ground  sufficiently  small  compared 
to  the  diameter  of  the  machine,  it  is  inherently  stable  in  attitude.  However, 
the  operating  height  required  for  stability  is  generally  incompatible  with 
the  height  needed  for  operation  over  typical  terrain.  The  range  of  stability 
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can  be  extended  by  using  interior  jets  which  effectively  "compartment" 
the  machine.  Since  the  individual  compartments  tend  to  exhibit  altitude 
stability,  the  machine  becomes  stable  in  pitch.  However,  as  shown  for 
example  in  Reference  5,  considerable  performance  loss  results  from  the  use 
of  internal  stabilizing  jets .  For  this  reason,  it  was  thought  worthwhile 
to  investigate  feedback  as  a  means  of  achieving  the  desired  degree  of  sta¬ 
bility  with  no  significant  loss  of  performance. 

The  dynamics  aspect  of  ground  effect  machines  has  been  considered  pre¬ 
viously  in  References  1,  2,  3  and  4,  although  not  primarily  from  the  point 
of  view  of  feedback  stabilization.  Reference  1  was  only  of  a  very  prelimi¬ 
nary  nature  although  feedback  control  was  considered.  The  final  results  of 
the  study  proposed  in  Reference  2  were  not  available  to  the  authors  at  the 
time  of  publication  of  this  report.  Reference  3  considered  only  machines 
which  had  been  previously  stabilized  by  means  of  internal  jets,  while  Refer¬ 
ence  k  considered  only  the  heave  motion. 

During  the  short  period  of  time  since  widespread  attention  was  first 
given  to  the  GEMs,  a  large  amount  of  work  has  been  done.  However,  a  number 
of  problems  still  remain.  This  is  particularly  true  with  regard  to  restor¬ 
ing  moments  and  damping  forces  and  moments,  since  most  of  the  attention  to 
date  has  been  confined  to  predicting,  and  experimentally  determining,  the 
performance  characteristics . 

A  short  review  in  three  sections  is  given  here  of  the  work  which  has 
been  done  . 

1.  Augmentation.'  A  study  of  the  general  proximity  effect  on  nozzles 
was  made  by  Von  Glahn  in  April  1957-  (Reference  6).  Shortly  thereafter  a 
theory  was  presented  by  Chaplin  (Reference  7)  which  assumed  the  jets  to  be 
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thin  and  neglected  the  effects  of  viscosity.  The  theory  was  extended  to 
thick,  viscous  jets  in  References  8  to  12.  A  qualitative  discussion  of  the 
effects  of  vortices  was  given  by  Nixon  and  Sweeney  (Reference  13).  A  large 
amount  of  experimental  data,  as  well  as  results  of'  experience  with  full 
scale  machines,  was  given  in  the  papers  presented  at  the  Princeton  Symposium 
(Reference  l4).  The  effects  of  vortices  were  treated  theoretically  in 
References  15  and  l6. 

2.  Moments.  A  theory  for  the  prediction  of  moments  employing  thin 

jets  was  given  in  Reference  17,  and  was  extended  to  thick  jets  in  Reference  18 . 
Several  experimental  investigations  of  the  static  moments  acting  on  GEMs  were 
reported  in  Reference  5* 

3.  Forward  Flight.  Experimental  investigations  of  GEMs  in  forward 
flight  were  given  in  References  19  and  20. 

Available  theoretical  and  experimental  data  seem  to  be  sufficient  to 
allow  prediction  of  the  augmentation  of  a  hovering  machine  with  good  accuracy. 
The  situation  with  regard  to  moments,  however,  is  not  as  good.  The  theory  in 
Reference  17  shows  the  machine  to  be  unstable  in  attitude  for  all  height  to 
diameter  ratios.  The  theory  of  Reference  18  gives  results  similar  to  those 
observed  experimentally.  However,  these  results  are  tied  to  certain  as¬ 
sumptions  and  the  sensitivity  to  these  assumptions  has  not  been  explored. 

In  summary,  actually  very  little  information  about  stability  derivatives 
is  available  in  the  literature .  Since  a  fairly  good  understanding  of  the 
vehicle  to  be  controlled  is  a  necessity  for  a  feedback  control  synthesis,  a 
large  part  of  our  effort  has  been  devoted  to  deriving  expressions  for  some 
important  stability  derivatives .  The  validity  of  the  theoretical  results  has 
been  checked  by  measurements,  and  the  results  are  used  in  a  preliminary  syn¬ 
thesis  of  a  feedback  control  system  for  stability  augmentation. 
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1.2 .  Equations  of  motion 

If  we  attempted  a  general  solution  to  the  problem  of  describing  the 
dynamics  of  the  ground  effect  machine,  a  system  of  six  coupled  nonlinear 
differential  equations  for  the  six  degrees  of  freedom  would  have  to  be 
solved.  However,  in  order  to  concentrate  on  the  dynamic  problems  which 
are  most  important  from  a  practical  engineering  viewpoint,  we  are  going 
to  make  simplifications  while  emphasizing  those  aspects  which  lead  to  a 
reasonable  formulation  of  the  control  problem.  Some  of  the  factors  which 
are  neglected  at  this  time  may  draw  more  attention  at  a  later  stage  of 
development  after  the  most  important  controllability  problems  have  been 
solved . 

The  first  in  a  series  of  simplifications  is  the  assumption  that  the 
longitudinal  and  the  lateral  degrees  of  freedom  can  be  considered  un¬ 
coupled.  The  most  obvious  coupling  terms  are  gyroscopic  moments.  A  dis¬ 
cussion  of  the  order  of  magnitude  and  the  scaling  of  the  gyroscopic  coupl¬ 
ing  is  presented  in  Appendix  A.  The  conclusion  is  that,  except  for  single 
propeller  ground  effect  machines  with  large  specific  loads,  the  gyroscopic 
coupling  is  practically  negligible;  therefore  it  will  be  neglected  in  the 
present  study.  Nevertheless,  it  should  be  mentioned  here  that  the  scaling 
of  the  gyroscopic  coupling  should  be  considered  if  the  dynamics  of  a  single 
propeller  ground  effect  machine  are  to  be  investigated  on  a  scaled  model. 

There  is  also  a  coupling  of  an  aerodynamic  nature  between  the  rolling 
and  pitching  motion,  due  to  the  fact  that  both  are  coupled  to  the  heaving 
motion  in  a  similar  way.  It  is  well  known  that  the  lift  does  not  remain 
the  same  when  the  GEM  is  not  parallel  to  the  ground;  the  lift  is  an  even 
function  of  the  attitude  angle  around  any  axis .  On  the  other  hand,  the  aero- 
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dynamic  moment  at  any  attitude  angle  is  also  a  function  of  height.  The 
mechanism  of  the  coupling  can  be  described  in  a  qualitative  way  as  follows. 

Assume  that  the  machine  is  at  a  certain  roll  and  pitch  attitude  angle  at  the 
same  time  and  that  the  roll  angle  is  changed.  This  causes  a  change  in  the 
lift  force,  therefore,  the  c.g.  of  the  machine  moves  up  or  down.  The  change 
in  height  induces  a  change  in  the  moment  around  the  pitch  axis  resulting  in 
a  pitching  motion.  Similarly,  a  change  in  pitch  angle  results  in  a  rolling 
mot  ion . 

The  strength  of  this  coupling  depends  mainly  upon  the  slopes  of  the  functions 
expressing  the  relationship  between  height  and  angle,  and  moment  and  height.  The 
effect  of  this  coupling  on  the  dynamics  is  strongly  influenced  by  the  damping 
in  heave  . 

The  even  function  L  =  $  (©)  where  ©  could  he  the  angle  around  any  axis, 
has  a  flat  maximum  at  0  =  0.  The  change  in  lift  with  small  deviations  from  the 
leveled  condition  has  been  shown  experimentally  to  be  very  small.  (See  reference  5 
and  Figure  3)*  Experiments  and  theory  (to  be  presented  in  subsequent  sections)  in¬ 
dicate  a  fairly  high  degree  of  damping  in  the  heave  motion.  Considering  also  that 
at  practical  heights  only  small  angles  are  possible,  this  coupling  between  the 
lateral  and  longitudinal  degrees  of  freedom  can  be  regarded  as  of  secondary  im¬ 
portance  and  is  neglected  in  the  subsequent  analysis. 

Other  coupling,  terms  may  appear  due  to  an  added  fin,  tail  or  auxiliary 
thrust .  Such  couplings  are  not  considered  here  because  we  are  concerned  pri¬ 
marily  with  low  translational  velocities  where  the  ground  cushion  exists  and 
where  the  aerodynamic  forces  exerted  on  a  fin  or  tail  are  not  very  strong. 

It  should  be  pointed  out  here  that  the  effect  of  these  additional  terms  on 
the  lateral-longitudinal  coupling  can  be,  and  should  be,  considered  carefully 
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in  the  case  of  any  particular  design.  We  intend  to  concentrate  on  the  dominant 
stability  problem  of  the  peripheral  jet  ground  effect  machine  at  and  near 
hovering . 

The  above  discussion,  as  well  as  experience  with  existing  full  scale  GEM 
models,  indicate  that  most  problems  of  stability  and  controllability  appear 
even  if  only  three  degrees  of  freedom  are  considered:  the  freedom  in  heave, 
in  attitude  around  one  axis,  and  in  motion  along  the  axis  perpendicular  to 
the  attitude  axis.  The  remainder  of  this  study  represents  an  approach  to  the 
analysis  of  the  stability  problems  for  the  purpose  of  a  preliminary  synthesis 
of  an  automatic  feedback  system  for  stability  augmentation. 

In  order  to  make  the  mathematical  treatment  of  our  problem  feasible,  we 
will  use  small  perturbation  theory.  This  approach  enables  us  to  investigate 
the  necessary,  but  not  sufficient,  condition  for  stability  of  a  nonlinear 
system;  the  system  must  be  stable  in  the  vicinity  of  any  equilibrium  point. 

This  investigation  reveals  approximate  quantitative  information  about  the 
most  important  stability  problems.  The  variation  of  the  stability  through¬ 
out  the  considered  flight  regimes  will  be  visualized  by  means  of  root  locus 
diagrams  so  that  basic  solutions  to  the  feedback  control  problem  can  be 
suggested . 

Stability  in  hovering  will  first  be  considered,  followed  by  a  discussion 
of  stability  in  forward  flight . 
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2 .  HOVERING  ANALYSIS 


2.1.  Hovering  equations 

As  an  approximation  to  the  stability  problem  in  hovering,  we  assume  that 
variations  in  attitude  are  sufficiently  small  so  that  any  lateral  or  longitudinal 
velocities  which  develop  are  quite  small,  and  therefore  forces  and  moments  due  to 
horizontal  velocities  can  be  neglected.  This  leaves  us  with  two  coupled  equations 
in  two  degrees  of  freedom.  Using  the  height,  V)  ,  and  attitude  angle  with  respect 
to  horizontal,  Q  as  perturbation  variables,  we  can  write  the  equations  in  linear¬ 
ized  form,  Laplace  transformed  with  zero  initial  conditions. 

For  heave: 


(-  Is  +  tr  s  ■*  e  ' 


£>L 

25© 


g>i- 

S5-, 


St 


and  for  pitch: 


(i) 


f£S.  <*> 

Sr  and  §&  represent  the  throttle  and  the  moment  control  input 
respectively. 

For  stability  investigations  only  the  homogeneous  equations  need  to  be 
considered.  After  dividing  by  the  weight,  W  ,  we  can  write  the  heave  equation 
in  the  following  form. 


s  - 


_L_  A  =  O 
Pl.  HQ  0 


(3) 


\-/S 


where 


The  moment  equation  can  be  normalized  by  dividing  by  L.D  where  D  can  be 
any  suitable  dimension  of  the  GEM 

(iss  aCM  h  -  O  w 

v  lp  ^  ae  z>o  /  ah 

In  effect,  this  procedure  serves  to  define  the  moment  coefficients. 

Let  us  consider  now  a  special  but  practically  very  important  case. 

For  equilibrium  at  the  leveled  condition  ^ =  0  since  Pu  is 
a  smooth  even  function  of  0  .  Also,  )0=o  =  0  since,  because 

of  symmetry,  the  moment  at  0  =  0  is  zero  independent  of  the  height.  There¬ 
fore,  in  hovering  with  zero  attitude  angle,  and  practically  also  at  small 
attitude  angles,  we  have  two  uncoupled  second  order  equations  of  motion  for 
the  two  degrees  of  freedom  in  heave  and  attitude.  This  enables  us  to  investi¬ 
gate  the  stabilities  of  the  uncoupled  modes  as  functions  of  height  and  this 
simplified  approach  will  lead  to  basic  suggestions  for  the  solution  of  the 
stability  problem  in  hovering.  The  effect  of  dynamic  coupling  and  forward 
flight,  will  be  discussed  subsequently. 

2.2.  Application  of  momentum  balance  to  GEMs . 

2.2.1.  Assumptions .  In  this  section  the  momentum  balance  approach  is 
used  to  derive  the  most  important  stability  derivatives  in  hovering.  The 
vortices  which  occur  because  of  viscous  effects  are,  to  a  certain  extent,  in¬ 
cluded  in  the  analysis  by  lumping  the  averaged  pressure  drop  under  the  base 
into  a  "base  efficiency  factor",  .  The  momentum  balance  approach  seems 

to  be  adequate  to  show  some  of  the  most  important  mechanisms  influencing  the 
stability.  The  relationships  to  be  obtained  are  also  very  useful  in  estimating 
the  effects  of  scaling  on  the  stability  derivatives. 
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The  momentum  balance  approach  has  been  used  previously  by  numerous  authors . 
The  method  given  here  differs  from  previous  presentations  in  that  basic  relation¬ 
ships  are  derived  using  only  the  momentum  balance  approach.  Actual  numerical 
results  may  differ  considerably  depending  upon  the  assumptions  which  are  made 
about  the  jet.  In  the  following  presentation  different  assumptions  about  the 
jet  may  be  taken  into  account  by  making  appropriate  substitutions. 

Some  fairly  general  relationships  concerning  the  base  pressure  and  the  air 
flow  at  the  periphery  of  the  base  can  be  deduced  from  the  assumption  that  the 
pressure  difference  at  the  periphery  is  kept  in  balance  by  the  sum  of  jet  moments. 
Following  Tulin  (Reference  25)  we  distinguish  three  different  regimes.  We  shall 
characterize  these  regimes  for  unit  length  of  the  periphery,  without  specifying 
any  method  of  determining  the  jet  momentum.  We  shall  assume  that  the  radius  of 
curvature  of  the  periphery  is  large  compared  to  the  height  of  the  machine  above 
the  ground,  i.e.,  we  assume  the  flow  at  each  point  of  the  periphery  to  be  two- 
dimensional.  (Three-dimensional  correction  will  be  made  later.)  We  shall  also 
assume  the  flow  to  be  incompressible  and  inviscid„  However,  viscous  losses  will 
be  considered  in  the  form  of  efficiency  factors. 

2.2.2.  Balanced  regime .  In  the  balanced  regime  the  jet  from  the  peripheral 
nozzle  seals  off  the  base  area.  At  the  peripheral  element  under  consideration, 
air  is  neither  entering  nor  leaving  the  base  area. 


For  horizontal  equilibrium  over  a  segment  of  the  periphery  of  unit  length. 


=  /^/S  +  fv*  ft) 


(5) 


where  K  =  the  height  of  the  nozzle  section  above  ground 

Tb  =  momentum  of  the  jet  leaving  the  nozzle 
3b.  =  final  momentum  of  the  jet 
=  a  jet  efficiency  factor 

If  the  base  of  the  machine  is  parallel  to  the  ground,  then  P-  Po 
where  |30  is  the  jet  inclination  angle  of  the  machine . 

Equation  5  can  be  normalized  by  dividing  by  the  jet  strength  out  of 
ground  effect.  Using  Bernoulli's  equation 

J  -  -V  v/*  -  *+  P_  <6> 


oO 


p  t  \^o  -  a  t  PT 


Dividing  (5)  by  (6)  we  get 


l 

a 


Pab  V.  _ 


=  *br( 


3b. 

T-o 


(7) 


Using  a  bar  to  designate  normalized  quantities  this  can  be  written  as:e 


5. 


P«b  P  =  /6  -t  "Jb.)  =  3b  (  \  t  /bun  j3  ) 


(8) 


where  the  length  is  normalized  by  the  nozzle  width  and  the  difference  be¬ 
tween  the  initial  and  final  jet  momentum  has  been  neglected. 


2.2.3.  Underfed  regime .  In  the  underfed  regime  air  is  flowing  out 
of  the  base  area  and  the  momentum  of  the  escaping  air  must  he  considered. 


Using  the  subscript  U-  for  the  underfed  regime  we  can  write  for  the 
horizontal  equilibrium 

■y  Pbu.  V>  =  h;r  (  3"u.  /^(3  +  Tuui-t  ^ 

=  VjT  ^_Tu.  O  *  £)  +  ^'^'2-1 

Here  Tu-i  is  actually  a  notation  for  the  normalized  change  of  momentum 
rather  than  for  the  jet  momentum. 

For  a  leve3.ed  machine  again  f3>=  (30  ;  if  the  base  is  at  an  angle  © 

to  the  horizontal,  |3>  =■  £>0  —  ©  .  The  mass  flow  of  the  escaping  air 

must  be  determined  from  an  additional  equilibrium  condition  for  the  entire 
base  area. 

2.2.4.  Overfed  regime.  In  the  overfed  regime  air  is  flowing  into  the 
base  area.  This  air  must  be  supplied  by  part  of  the  nozzle  jet.  This 
weakens  the  total  jet  momentum  due  to  diminishing  of  the  mass  flow  which 
is  turned  away  from  the  base  area  and  due  to  the  momentum  of  the  air  flow¬ 
ing  into  the  base  area. 


18 
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For  horizontal  equilibrium 

-j  Pso'v,  -  Ht  (To  £  +  Toi  -  To*) 


(10) 


where  again  3  =  3o-  ©  .  The  mass  flow  of  the  air  flowing  into  the 

base  area  must  again  be  determined  by  means  of  an  additional  equilibrium 
condition  for  the  base  area. 

Using  the  basic  concept  of  these  three  regimes,  equations  for  forces 
and  moments  can  be  determined  in  a  form  still  independent  of  the  particular 
choice  of  jet  theory,  but  nevertheless  indicating  the  basic  relationships 
and  important  factors .  It  should  be  noted  that  the  basic  concept  can  be 
used  also  for  other  than  simple  peripheral  jets. 

2.3.  Lift  and  augmentation. 

The  over-all  lifting  force  acting  on  a  ground  effect  machine  with  a 
single  peripheral  jet  can  be  expressed  in  the  following  general  form: 

L=  3  cm  1  J  PB  (11) 
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where  3  is  "the  jet  momentum  per  unit  length  of  the  periphery,  jL  is  the 
length  of  the  periphery  and  $  is  the  base  area  enclosed  by  the  jet .  In 
a  leveled,  balanced  condition  the  height  of  the  jet  is  uniform  along  the 
periphery  and  0=  0O  ;  the  base  pressure  is  considered  uniform  over 
the  whole  base  area,  L  is  an  efficiency  coefficient  expressing  the  ex¬ 
perimental  fact  that  due  to  vortices  and  other  losses  only  part  of  the 
computed  pressure  contributes  to  the  lift.  For  the  leveled  balanced 
condition 

L-  rlT  Jb  l  ow-iio  t<lLSPe  (12) 


Generally  both  Tb  and  Ps  are  functions  of  the  height,  V~)  .  We 

introduce  now  the  specific  load  Pu=  L/S 

Pl_On)=-  Fjj  CV>)  -4  C04.  0o  -t  Hl.  (^) 


For  determining  the  augmentation  factor  the  reference  is  the  specific 
lift  produced  by  the  same  total  head  out  of  ground  effect-  with  0O  —  O 
and  the  same  jet  efficiency  factor,  .  Out  of  ground  effect  Fg =  ° 

and  with  Equation  6  the  augmentation  can  be  expressed  as 


A(v >)= 


a  nT 


Pu  IF)  _s 
Pr  It  “ 


=  3b(F) 

3"  oo 


Ol 


_S_ 


P^Ch) 

Pt 


(i^) 
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The  total  derivative  is  used  from  here  on  since  we  are  considering  freedom 
in  height  only. 

The  specific  restoring  force,which  is  in  effect  the  normalized  "spring- 
constant"  in  heave, is  defined  as 


ft* 


-  Aft- 1 
dL  Vi  '  h  = 


(19) 


cLTb  JL 

dLb  S 


0*4-  $o  +  Hl 


cLPb  1 

cl  W  -1  h  =  he 


The  derivation  below  is  based  on  the  following  physical  considerations. 

1.  Among  the  related  variables  V,  ,  p6  ,  Jb  there  is  only  one 
independent  variable . 

2.  With  our  assumptions  the  jet  is  determined  uniquely  by  the  normal¬ 
ized  base  pressure,  Pg  ,  therefore 

cL  Tb  _  A.  3  b  cL  Pb  (20) 

d.h  cLP6  d-b 


The  change  in  base  pressure  with  height,  P_  ,  can  be  determined  from 

AM 

Equation  5  hy  differentiating  and  substituting  the  above  expression  (20). 
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Substituting  Equation  20  and  the  above  expression  into  Equation  19,  we  get 


k-  *1u  t: 


i  •+  'Jz 

!l  Wn . 


Dl  d  Jb  n 

WZ  dTPa  8 


I  - 


if  ( *+/**&.) 


and 


normalizing  by  T<o  =  A  t.  Bp  with  Cn  = 


K  * 


Pb 

n-  17 


1  *  Se 


(22) 


=  n. 


fj?  ^b 
Vi  <=h 


This  equation  also  serves  as  a  definition  for  C  *>  and  C 

The  fraction  in  which  appears  both  in  the  numerator  and  in 

d  Pe 

the  denominator  indicates  the  influence  of  the  change  of  jet 

momentum  with  base  pressure .  If  the  jet  would  not  change  with  the  base 
pressure,  f"g  ,  the  spring  constant  would  be  simply, 

*  -  Hu  \  (23) 
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jfhe  correction  term  in  the  denominator  expresses  the  fact  that  this  full 
cannot  develop.  For  example,  downward  motion  from  equilibrium  results  in 
some  weakening  of  the  jet  and  less  restoring  pressure  can  develop  than  with 
an  unchanged  jet.  The  correction  term  in  the  numerator  expresses  the  change 
in  the  direct  contribution  of  the  jet  to  the  lift . 

2.4.2.  Damping  in  heave .  If  the  machine  is  moving  downward  with  a 

•  ^ 

velocity,  h  ,  air  must  leave  the  base  area  and  if  the  machine  is  moving 

with  an  upward  velocity,  air  must  flow  into  the  base  area.  Damping  forces 
arise  because  of  changes  in  the  over-all  lift  at  any  height  which  occur  in 
either  an  underfed  or  an  overfed  regime.  The  mass  flow  leaving  or  entering 
the  base  area  is  determined  by 


m  JL  = 


Sh 


(24) 


where  IT)  is  now  the  mass  flow. 

For  simplicity  only  the  leveled  condition  is  considered.  Here  the  jet 
is  uniform  along  the  whole  periphery. 

2. 4. 2.1.  Downward  velocity.  While  losing  altitude  with  a  velocity 

O 

V~>  4  O  ,  the  escaping  air  forms  an  additional  jet  with  momentum  Ju.1 
over  a  unit  length  of  the  periphery,  resulting  in  an  underfed  regime. 
Because  of  the  leveled  condition,  P>-  .  Indicating  the  changes 

from  the  balanced  condition  we  obtain  from  Equation  9 

— —  ^  ^Bu.)  Vi  =  +  ■+  3  (^5  ) 
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We  get  the  pressure  difference  due  to  the  escaping  air  if  we  subtract 


Equation  8  from  Equation  25 . 


1 


A  Pgu.  V'l  .  =■ 


A  Tu,  (  \  -t  /oon  (3c)  + 


0^2. 


(26) 


_  cL  Tb 

T5 


8 


A  Peui.  C  1  (3o)  +  "3  u-i. 


The  predominant  term  on  the  right  hand  side,  JianL  ,  indicates  a 
pressure  increase .  Alu  arises  because  the  mass  flow  at  the  nozzle 
is  a  function  of  the  base  pressure.  An  increase  in  base  pressure  causes 
a  decrease  in  mass  flow,  therefore  ATTul.  is  a  negative  number.  Solving 
Equation  26  for  APga  we  obtain. 


A  P eu.  = 


Ot 


\  -  A  \  +  /U*\  Pc') 


h 


(27) 


^  jQ  T  Tu.2, 

C.  V-i  Vi 


where  Gy-,  is  the  same  as  defined  in  Equation  22. 

Differentiating  the  specific  lift,  Equation  l6,  and  using  the  base 
pressure  change  expressed  by  Equation  27,  we  obtain  the  following  ex¬ 
pression  for  the  normalized  damping  force, 
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APL 


u 


=  (A  nT  O-  co4./6o  +-  r\L_)  aP8 


(28) 


•  -  Ani_nT  % 

t-b  n 

Cb 

where  -=-  has  been  defined  in  Equation  22 . 

L  b 

The  additional  momentum,  Ju.a  ,  can  be  determined  from  the 
mass  flow,  Equation  2k,  and  Bernouilli’s  equation 


Ju/2=--£  \^/>pe  h  Vi  o 


(29) 


In  normalized  form 

S  I 


Tu,z  =  -  x  ^pt  v,  =  -  -k  \TS~ 


o  Jo. 


(30) 


where  is  the  jet  velocity  out  of  ground  effect. 


Substituting  into  Equation  28 


aPl 


dL*n 


£  ^  Ht  Cb  4-  xTpe  — 

o-  Cb  h.  ’  ^r- 


(31) 


2A.2.2.  Upward  velocity.  While  gaining  altitude  with  a  velocity 
h  >  o  j  air  is  being  supplied  by  the  jet  into  the  base  area,  resulting 
in  an  overfed  regime.  For  the  leveled  machine,  /*=  (So  .  Indicating  the 
changes  from  the  balanced  condition.  Equation  10  can  be  rewritten  as  follows 
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A  Pbo)  V)  =■  IQj  jj^Tb  A  To)  /Ujo  (3o  -»■  Toi  -  ToiJ  (32) 


We  get  the  pressure  difference  due  to  the  flow  into  the  base  area  if  we 
subtract  Equation  8  from  Equation  32: 


A  Peo 


=  'Or  \^A.To  £>o  -(^Tbi  -  Toi)  -  To2^ 


We  can  arrive  at  a  simple  approximation  to  the  right  hand  side  as  follows: 


Tbi  =  m0  ('To0  where 


> 

\T~ 


Tq2  =  ft)  Woi  where  <VJoa  =  V  FS~  -  Pe 


where 


J01-  (<™0  +  Am0-m)  (VToo  A  moKToo=  aTqi  =  aT< 


If  we  normalize  by  Too =  HOoo  Woo  where  m co  is  "the  mass  flow  at  the 
nozzle  out  of  ground  effect  we  get 


Tbt  ~  To  1  — 


PQ  _  A  -lot 
Too 


=  m  -  ^Tc 


I02  =  'no 


\\  - 
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(§> 


Using  again  aT=  4*  ■A  Pa  >  ve  find,  by  substituting  into  Equation  33 


dLPs 


APeo  =  -  m  (  I  +  J  i-  PB  ) 


(35) 


where  G^  is  the  term  as  defined  in  Equation  27,  and  using  the  mass 
flow  Equation  24 

(36) 


APBO  =  -  ^  ^ 


h_ 

VoO 


V>>  o 


The  change  in  the  specific  lift  describes  the  normalized  damping 


A^=  cbnu^p6o=  Is 


(37) 

K3V>o 


For  determining  a  damping  factor  for  comparison  with  experimental  data  we 
consider  the  average  damping  force  acting  over  a  complete  cycle  of  oscillation. 


I"  + 


rk^Cb. ( i  +-  +  J » —  Pe ) 

CX  h  ew  v  / 


IhJ 

UoO 


(38) 


We  make  an  error  of  less  than  4$  for  Q«l  <  PB  <C0.S  ,  i.e.,  for  all 
practical  values  of  P6  ,  if  we  write 
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j-  DaP«—I+  lAp-«l]=  TT  t 


(39) 


This  is  an  expression  for  the  normalized  average  damping.  We  obtain  the 
actual  damping  in  terms  of  pressure  if  we  multiply  both  sides  by  Ry 


t 


(4o) 


-  \.LS  jIS  \V)\ 

h 


ex. 


2.^.3.  Frequency  and  damping  ratio  of  the  heaving  motion.  The  transformed 
differential  equation  describing  the  heaving  motion  with  the  coupling  term 
neglected  can  be  written  as  (See  Equation  3): 


c 


—  2.  4r  ^r.  O 


_  -St  dR_ 

Pi_  Ah 


)*- 


(4l) 


Using  Equations  19  and  22,  the  undamped  natural  frequency  is, 


OO 


H-  If 


h  'J 


Hi,  Pgg  db 
Pu  Ch 


(42) 


a/t 


V>  'J 


Rl._Pr  eb 

Pc.  Ch 
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We  can  determine  the  relative  damping  ratio,  ,  from  A  5"h<^>h=  ^ 

Pi_  Oh 

Expressing  as  the  coefficient  in  Equation  40 


0.f?-5 

CX, 


(43) 


The  influence  of  scaling  is  an  increase  of  the  damping  ratio  with  the 
square  root  of  the  linear  dimension.  The  volume  of  air  which  moves  into 
or  out  of  the  base  area  increases  proportionally  with  the  base  area,  whereas 
the  "nozzle"  area  through  which  this  air  moves  increases  proportionally  with 
the  length  of  the  periphery. 

The  damping  ratio  also  varies  inversely  with  the  square  root  of  the 
specific  load,  P u,  .  The  total  pressure,  Pj-  ,  must  increase  with  P,_ 
and  the  base  pressure,  P5  ,  is  proportional  to  Pp  •  Therefore  the 
variation  of  the  velocity  of  the  air  moving  into  or  out  of  the  base  area 
is  reflected  in  the  factor. 

2.5.  Pitching  motion. 

2.5.1.  Moments .  The  physical  basis  for  the  derivation  of  an  expression 
for  the  moments  acting  on  a  ground  effect  machine  as  given  below,  is  the  same 
as  the  one  adopted  by  Lin  (Reference  18)  and  Webster  (Reference  IT ) •  The 
following  derivation  is  different,  however,  in  the  respect  that  a  more  general 
result  is  obtained  so  that  the  effects  of  different  mechanisms  and  various 
assumptions  can  be  discussed. 

Consider  that  the  machine  is  in  equilibrium  at  some  small  attitude  angle. 
The  experimentally  observed  flow  picture  is  similar  to  the  one  shown  in  the 
following  Figure . 
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The  cross  flow  is  observed  experimentally  and  can  be  explained  by  con¬ 
sidering  the  balance  of  pressure  and  momentum  forces.  When  the  machine 
is  tilted  to  the  right,  for  instance,  the  force  due  to  the  base  pressure 
has  a  component  parallel  to  the  ground.  The  machine  must  thus  exert  an 
equal  and  opposite  force  on  the  fluid  mass  under  the  base.  In  order  for 
this  fluid  mass  to  remain  in  equilibrium,  the  pressure  force  must  be 
balanced  by  a  momentum  change  plus  the  frictional  force  at  the  ground. 

Since  the  skin  friction  force  is  quite  small,  the  pressure  must  be  bal¬ 
anced  mainly  by  a  momentum  change  which  can  only  be  brought  about  by  part 
of  the  air  from  the  low  jet  flowing  toward  the  high  jet. 

The  losses  due  to  viscous  effects  including  the  pressure  drops  caused 
by  the  vortices  are  lumped  into  the  lift  efficiency  factor,  .  Diff¬ 

erential  pressure  changes  caused  by  the  difference  between  the  vortices  at 
the  two  ends  are  considered.  The  result  of  the  derivation  will  show  that 
these  assumptions  are  sufficient  to  explain  the  nature  of  the  moment  derivative. 
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and  the  results  will  be  later  modified  for  application  to  a  three 
dimensional  machine. 


Let  us  take  a  section  of  the  vehicle  of  unit  width  and  denote  the 
distance  between  the  two  jet  sections  by  ^  f'  .  If  the  base  is  inclined 
with  respect  to  the  ground  at  an  angle,©  ,  as  illustrated  in  the  above 
Figure,  there  is  an  underfed  regime  at  the  high  end  and  an  overfed  regime 
at  the  low  end.  The  mass  flows  into  and  out  of  the  base  area  must  be.  equal. 
Our  unknowns  are  aP0l=  P@u -  PBe.  ,  A  Pen  =  Pew  -  Fee-, 
and  the  mass  flow,m  .  We  have  two  momentum  balance  equations  for  the  two 
ends;  the  third  equation  is  provided  by  writing  the  horizontal  momentum 
balance  equation  for  the  total  base  area.  Linear  pressure  distribution 
along  the  base  will  be  assumed. 

For  the  high  end  we  can  write,  using  Equation  5,  and  the  notations  of 
the  above  Figure . 

-  -  P  -  -r  (WO 

p8H  =  nT  |_  JH  (.  '  +-  j 3  w)  •+  Tui  j 
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(45) 


_J_  (pee ^  )(  he  +  r>  ©^)  ”  Ht  aTh)  D  't'^MA  (Po  ©)I 

+  XhA  ■+  ^”fc>  +  -A  Th  ,} 

We  now  make  small  angle  approximations  for  ©  and  neglect  the  second 
order  terms 

Pge  Vtg  -+  Pbs  r6+  W  aPbh  +-  P  ©  A  P  BH  - 


=  ^Hn-  [(Tb +aTh)(  '  •V ^AJUV\  |3o)  —  ■+  ATh)&  Otr<L  |3o  +  Th*"] 

Subtracting  the  equilibrium  equation  for  the  leveled,  balanced  condition, 
p bcYyz  =  2.  Y\t  Tto  (  \ -t/AM  (3o)  and  substituting  aTh'  A  PBh 

we  get 

C>e+  C'0)aPBh  =  (4T) 

=  “  PW  r©  2-  "H-x  ^  (.'+■  aPbw  -  Xb©  OJr^-  Po-*"  ^ 

where  P'  =  P  +  ^p&  C'<"L 

Rearranging  and  using  Ch  as  defined  in  Equation  22  the  differential 
momentum  balance  for  the  high  end  can  be  expressed  as 

(GhW  +Ti©)aPbw  =  (,+8) 

=  -  Pbci"P  ©  -  >  ~^T  Tb  ©  O&o-  +  <2  r^j-  m  (  Nr  Ha-  ^ru-2.') 
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where  is  the  normalized  cross  flow  and  cK  (0"|_q_  is  the  normalized 

velocity  of  the  cross  flow  at  the  high  end. 

The  left  hand  side  of  this  equation  expresses  the  pressure  changes. 

The  three  terms  on  the  right  hand  side  indicate  clearly  the  three  major 
factors  contributing  to  the  pressure  change.  The  first  term  is  due  to  the 
increase  of  the  nozzle-to-ground  distance,  the  jet  momentum  can  balance 
only  a  smaller  pressure  over  a  larger  area.  The  second  term  is  due  to 
the  change  in  jet  incidence  angle  at  the  nozzle  with  respect  to  the  ground, 
decreasing  the  horizontal  component  of  the  jet.  The  first  two  terms  are 
of  stabilizing  nature,  the  third  term  which  is  due  to  the  cross  flow  is  a 
destabilizing  term. 

For  the  low  end  we  can  write,  using  Equation  10,  and  the  notation  of 
the  Figure  on  page  32 

Pbl.  hu  *  z3*- 

or 

(FW  +•  he 

=  Htt  a T A0’0  C  ©)  3ui  _  Tu-z  1 

We  perform  a  sequence  of  steps  similar  to  the  manipulation  of  the  equation 
for  the  high  end,  combined  with  the  steps  taken  previously  to  express  the 


conditions  for  the  overfed  regime. 


(50) 
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and  making  similar  steps  as  before 


—  "ili  [Ub-*  a3h)  (/ua\|30  -  ©  c-oo.  (3^)  -(Tb-t  aTu)(a^»|5c.  +  ©c«Ae) 
■+  3Vu  ■+  ThJ.- 

=  P|j  L-a  Jb  ©  ci^f3o  -t  ^5*  CaTP0m-  aPbl_)  /uvnySo 

^Jb  0  o  (a  Pen  +  Aw  Psu\  •+  Tb  t  A3h 

dL  Ps  ' 

Th2—  Ti_|  1 

Substituting  —  3»_|  «  Ch  -  ATl.  ;  VT>  NTVl2.  J 

and  3,  (,  V  i-  /U/Vn  =  P>e  (  l  —  G  b) 

We  arrive  at  the  following  final  form  of  our  third  equation 

Cnur'0  -  (  » -Gh)helApBH  +  [riu^'©  +  (»  -ChJhe^AFW  (53) 

=  Pee  r  ©  -  4-  i1TTb  ©  gcxl  NrH2^) 
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Since  the  moment  will  be  proportional  to  A  Pgv-t  —  A  Pbl_ 

it  is  advantageous  for  the  calculation  to  choose  this  as  one  of  the 


variables.  Multiplying  Equations  48  and  51  by  and  writing  their 


difference  and  their  sum  we  arrive  at  the  following  set  of  equations. 


V\l  V"'©  (  APen  ■+  a.Pbl)  -  (  1  -Ch)  C  APeM“  APeL) 

—  2.Y)T  (_  1  -v  (VTm  2.^  ^  ~  <3  (o  t_  Pse  r  -V  A  Ht  ~^t>  ccd^o^© 


Hu^hVie  (a  Psh  4  APg^  +  Y" |i_r'©  (aPbh-  aPbl) 

\  -t  i?u  -  o(  Ntl.z)  -  O 


(54) 


Hl.  'r<  ©  APbl)  Y'lu  Qh  ^>2.  (a  Pew  —  APcl) 

— Hi.  Plj  rr\  (  '  ^hz  —  <*KTu2.^~ 

=  -  A  (nL  psc  r  -t  5.  riuViT  3b  ^|9o)  © 


3T 


The  solution  for  A  Pen-  A  PBU  can  be  expressed  in  the  following 
form: 


A  Pew  -  A  P©i_ 


Q.(  Q  q2 
c3  -V  +  C.S  ©2 


(55) 


We  wish  to  determine  the  moment  derivative 


dLKl 

dL© 


at  ©=  o  . 


This  restriction  is  justified  because  experiments  indicate  that 


dLM 

a© 


is  fairly  constant  in  the  vicinity  of  zero  degree.  Therefore  we  need  to 
determine  only 

(56) 


a© 


(aPsh  -  A  Peu)  \  »  -§i- 


0  =  0  C3 


he  £ 


r  -v  a  Ht  Tb  cj&4.  ^0)  tou.  (v-ot)  -  (\  — h|i_)  Pqo  r  ( v  r*a)~] 

I  -+  Whi  +  C  1-  Q*)  \TU^  C  t  -c*  )3 


The  moment  acting  on  a  segment  of  unit  width  is  easily  derived  as  follows 
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Assuming  linear  pressure  distribution. 


Pr  — 


P&H 


~  ^  PgL-^) 


aM  =  zipb  *.  dx 


Integrating  this  gives*, 


A[  Ap8^^  =  J  CaPQH-APbl)^’ 


We  normalize  by  ^  L.  f1  =  4-  Pt_  r2  where  L.  is  the  lift  acting  on  the 
segment  and  the  area  is  S  —  f'  .  The  moment  coefficient  of  our  two 
dimensional  segment  is 


at-r  12  R_  (A Pbh“ aPBu)=  (aPBv_, 


^Pbm- aP0u) 


The  moment  derivative  at  ©=  o  is 

<ACm  =  <  <L  r aPbh -aPBl)  = 

d©  l2Pi_ 

_L  £  (  nu  Pee -t- a  T- CM-fto")  NT-ut.  (\-oQ-(>-*W)  Pee-  0-^*0 

^Pu  *"*  '  **  (,  l  -  Cva)  KTu'l-  C  '  ~ ) 


*There  is  an  additional  moment,  the  change  in  jet  momenta  multiplied 
by  their  moment  arms .  It  can  be  shown  that  this  moment  is  negligible 
as  long  as  r-  / 1  »  \ 
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and  that 


Recognizing  that  in  Pu  =  £  Tb  1^|T  o_  C^r4^  ■+  Hu 

at  i 

for  our  segment  Cl=  —  «  -=-  ,  we  can  write  the  moment  derivative  in 

the  following  form: 


=  (60) 

d© 


x. 

< o  he 


\Tl2.  C  I  -  oQ 

I  -t  Kt  h~2. 


-  (  I-Hl) 


P&Z 

PL 


'♦i 


dLTb 


(vrur  (h-«)  (  ^ 

4  WTK2  ^ 


/3o^) 


This  result  clearly  indicates  some  of  the  influences  of  different 
physical  parameters  and  assumptions  on  the  moment  derivative.  The  two  terms 
in  the  numerator  of  the  coefficient  of  r/a*  are  of  distinctly  different 

origins.  The  positive  term  1  ~  °*  ^  could  be  called  a  cross  flow  variable. 

If  WHi 

The  parameter  KX lz  is  the  normalized  velocity  of  the  inward  flowing  jet 

at  the  low  end,  and  it  can  be  approximated  by  ATl.2.=  \  V  ~  Pee  since 

this  jet  is  flowing  into  an  area  where  the  average  pressure  is  Pb  .  As 
to  fcr M2.  the  picture  is  not  so  clear.  We  know  that  the  mass  flow  must 
be  the  same  as  that  of  the  inward  flow  at  the  low  end  and,  therefore,  the 
change  in  jet  momentum  at  the  high  end  is  determined  by  the  change  in  velocity 
of  the  outflowing  jet.  Because  of  dissipation,  the  cross  flow  will  not  have 
the  velocity  U2.  when  it  reaches  the  high  side  .  We  designate  the 

velocity  of  the  cross  flow  when  it  reaches  the  high  side  by  o  ItT u"2- 


where  C<  is  a  parameter  which  can  have  any  value  between  zero  and  one . 

Using  BernouiJ.li 's  equation,  we  than  find  that  the  normalized  velocity 
of  the  cross  flow  outside  the  base  area  is : 

AJna  =  jPe  +  a1  (  I-  Pe) 

When  there  is  complete  dissipation,  CA  —  O  and  we  obtain  the 

minimum  estimate  for  to"  Ha 

to”  H2  = 

If  there  were  no  dissipation  CX  would  be  equal  to  one  and  this  would 
correspond  to  expansion  from  the  total  pressure  to  the  ambient  pressure  and 
would  give 

-  I 

The  quantitative  effect  of  different  values  of  will  be  shown  later. 

The  positive  "cross  flow  term"  is  destabilizing.  If  the  other  numerator 
term  were  zero,  this  term  alone  would  make  the  ground  effect  machine  statically 
unstable  at  all  heights. 

The  second  and  stabilizing  term  could  be  called  a  "base  efficiency  term." 
The  (  \  —  )  factor  indicates  that  this  term  would  vanish  if  there  were 

no  losses  under  the  base.  We  have  incorporated  in  all  losses  resulting 

in  a  pressure  loss,  including  the  averaged  pressure  drop  due  to  the  vortices. 
Therefore,  we  can  state  that,  in  an  indirect  way,  the  vortices  provide  the 
source  of  a  stabilizing  mechanism. 


We  have  not  considered  any  differential  effects  of  the  vortices;  hut 
such  effects  can  be  expected  to  be  of  secondary  importance  if  our  assumptions 


lead  to  a  good  quantitative  agreement  with  experiments.  That  this  is  the  case 


will  be  shown  later. 

As  to  the  denominator  the  consistently  negative  sign  of  the 


makes  the  sum  in  the  denominator  always  less  than  one , 


.  ..  dL  T0 

of  the  - 

A.  PB 
L_~2  l  \  -~d 

\  ■+■  kr hi 


factor 


decreases  as  Ve  decreases,  so  that  their  ratio  is  not  expected  to  vary 


by  any  large  amount .  However,  at  very  low  heights  ^ 

dL  r© 


can  be  expected  to 


increase,  thereby  causing  a  decrease  of  the  denominator  and  an  increase  in  the 
absolute  value  of  the  moment  derivative.  The  physical  meaning  of  the 
factor  is  that  it  expresses  the  contribution  of  the  changes  in  jet  momentum  to 
the  moment . 

It  is  interesting  to  notice  that  both  the  stabilizing  and  the  destabilizing 


terms  in  the  moment  derivative  have  their  physical  origin  in  the  losses  occurr¬ 


ing  in  the  base  area,  according  to  the  momentum  theory.  If  the  lift  efficiency 
factor,  V\\_  ,  were  unity  there  would  be  no  stabilizing  moment.  If  the 

cross  flow  were  not  dissipated  at  all  (  CX  =  i  ),  there  would  be  no  destabilizing 
moment.  Consequently,  if  no  losses  occurred  in  the  base  area  the  moment  at  0=0 


as  well  as  the  moment  derivative  would  be  zero  at  any  height. 

2.5.2.  Three  dimensional  corrections.  We  now  have  to  make  an  estimate  of 
the  correction  to  be  used  for  the  application  of  our  moment  derivative  result 
to  three  dimensional  machines. 

We  shall  normalize  the  moment  of  the  three  dimensional  machine  by  dividing 
the  total  amount  by  where  L.  is  the  total  lift,  S  is 

the  total  base  area  and  is  the  maximum  half-length  (half -width )  of  the 

machine .  As  an  approximation  in  the  following  derivation  we  shall  neglect 


the  Tb  term  in  the  moment  derivative .  This  term,  which  represents  the 
direct  contribution  of  the  jet  momentum  to  the  total  moment,  is  small  at 
practical  heights.  Using  the  notation  of  the  Figure  shown  below  and  re¬ 
ferring  to  the  derivation  of  Equation  59>  we  can  write 


fcPi_SRz 


I  +  nr ^2. 


-  (  V- n«_)  Pee] 


1  <L Jw  ,  \  nr.-iii-a) 

' +  -frszi 


r-J  dr, 
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Comparing  this  result  with  Equation  59>  the  first  bracket  can  be  recognized 
as  the  two  dimensional  moment  derivative,  except  for  the  omitted  3"t»  term. 
The  second  bracket  is  the  correction  factor.  Applying  this  result  to  the 
entire  moment  derivative,  including  the  3"b  term,  we  can  write 


A.  Cm  \ 

d-G  1 3D 

where 


r 


cLCjm 

&e  2P 


(62) 


This  correction  factor  is  essentially  the  same  as  that  obtained  by  Webster 
(Reference  17 )  who  made  somewhat  different  assumptions.  The  factor 
is  the  ratio  of  the  rectangle  enclosing  the  periphery  to  the  base  area. 

The  integral  represents  the  effect  of  the  shape  of  the  machine  on  the 

correction  factor.  For  a  rectangular  machine  r-  <  ,  for  a  circular  machine  r=  .ns 
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It  should  be  considered  that  there  is  no  practical  need  for  a  greater  accuracy 
in  estimating  the  moments.  Therefore,  our  approximation  to  the  correction  factor 
by  simply  using  the  area  ratio  can  be  considered  adequate  for  our  purposes. 

2.6.  Pete rmi nation  of  base  pressure  and  jet  momentum. 

In  all  of  the  expressions  given  previously  for  frequency  and  damping  of 
the  heave  motion  and  for  the  pitching  (or  rolling)  moment,  it  is  necessary  to 
know  the  base  pressure,  Pg  ,  and. the  jet  momentum,  Tgj  ,  as  a  function  of 
the  height,  Pi  ,  and  the  physical  parameters  of  the  machine.  A  number  of 
different  expressions  can  be  obtained  for  Pg  and  Jg  depending  on  the 
assumptions  made  about  the  flow  pattern  or  the  jet  momentum.  We  will  briefly 
review  three  methods  of  obtaining  these  expressions  which  have  been  given  in 
the  literature . 

1.  Thin  jet  theory  (Chaplin,  Reference  7) 

The  balance  between  pressure  forces  and  momentum  change  for  a  level  GEM 
can  be  expressed  as, 


PbOA)«  j  +  (3a)  (  63  ) 

where  jl  is  the  length  of  the  jet.  If  the  jet  thickness,  "t  ,  is  small 
compared  to  the  height,  \n  ,  we  can  assume  the  momentum,  T  ,  is  constant 
and  equal  to  its  value  out  of  ground  effect.  Then, 


T=  AJL-t  (  '/a  p  Vj")  =  aUPt 

where  V*  is  the  velocity  of  the  jet.  Therefore, 

Pb  ~  ^  =  a-(  l  =  ^ 

which  defines 


(64) 


(65) 
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Obviously,  this  theory  cannot  be  expected  to  yield  good  results  near  and 
within  the  height  range  where  it  gives  Pg  >  |  and  this  range  is  within 
Gour  range  of  interest . 

2.  Exponential  theory.  (Stanton -Jones,  Reference  12) 

In  this  approach  we  assume  that  Equation  63  is  satisfied  for  a  diff¬ 
erential  element  of  the  jet,  i.e., 

P>dP  =  T  (  l  -V-  (lur\  (3)  (66) 

where 

t  =  p  at 

and  V  is  the  jet  velocity  at  the  element  under  consideration. 

From  Bernoulli! 's  equation, 

_i_pv2=  pt-p 

Combining  these  equations  and  integrating  across  the  jet  gives, 

P©  -  \  -  or**- 


(67) 


(68) 


Substituting  the  Taylor  series  for  the  exponential,  we  note  that  this 
approaches  the  value  given  by  Equation  65  for  small  fyC  . 


Pe  -  /y. 


CmO3 

a  1  3  1 


b6 


3.  Uniform  jet  theory.  (Fuller,  Reference  9;  and  Rethorst,  Reference  11) 
If  we  assume  that  the  jet  near  the  machine  has  approximately  constant 


thickness,  then  the  balance  of  pressure  and  centrifugal  force  in  the  jet 
gives 

d£  =  Q_s/*  (69) 

dp  /  p 


Using  Bernoulli's  equation  and  integrating  gives, 


P 


8 


t  - 


(70) 


where  Po  is  the  radius  of  curvature  of  the  outside  of  the  jet  and  P, 
that  of  the  inside.  Expressing  these  in  terms  of  height,  Vl  ,  jet  thick¬ 
ness,  ”t  ,  and  jet  inclination  angle  (5  0  and  defining 


\n  —  1 

V\  •+■  i&jf1  f3  o 


gives, 

Pe  =  l  ~ 


(71) 


In  the  thin  jet  theory  it  was  assumed  that  the  jet  momentum,  , 

constant.  In  the  last  two  theories  we  can  easily  obtain  the  following 
expression  for  *3"t>  (i.e.,  divided  by  3" y-,  out  of  ground  effect) 


~r  _ \  (  -  2.  (V-  N 

■^b  ~  V  '  ^  ) 


(72) 


k7 


was 
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3.  DISCUSSION  OF  THE  THEORETICAL  AND  EXPERIMENTAL  RESULTS 


In  the  preceding  section  theoretical  expressions  were  obtained  for 
several  of  the  important  stability  derivatives  for  the  hovering  GEM. 

The  variation  of  these  derivatives  with  changes  in  the  important  physical 
parameters  are  shown  by  means  of  a  number  of  plots  given  at  the  end  of 
the  report.  These  results  are  discussed  in  the  following  section  along 
with  results  obtained  from  a  series  of  experiments  with  an  eight  foot 
diameter  GEM  model.  The  experimental  set  up  and  the  procedure  used  for 
obtaining  the  data  is  discussed  in  detail  in  Appendix  B. 

It  was  shown  previously  that,  for  a  vehicle  initially  trimmed  paral- 
el  to  the  ground,  the  heave  and  pitch  motion  can  be  considered  uncoupled. 
We  thus  divide  the  discussion  into  two  sections:  heaving  motion  and 
pitching  motion. 

3*1.  Heaving  motion. 

We  previously  obtained  expressions  for  the  variation  of  the  lift 
force  with  height  and  rate  of  change  of  height  and  used  these  results  to 
find  expressions  for  the  frequency  and  damping  ratio  of  the  heaving 
motion.  Figures  i-  and  5  are  plots  of  the  undamped  natural  frequency  of 
the  heave  mot'.c  times  the  square  root  of  the  jet  thickness  (in  feet) 

COh  {V  (see  Equation  b2) ,  as  functions  of  jet  inclination,  (3 o 
nozzle  to  base  area  ratio,  Ou  ;  lift  efficiency  factor,  L  .  The 
general  character  of  these  curves  is  similar  to  the  curve  obtained  by 
thin  jet  theory  except  for  low  heights.  In  the  thin  jet  theory,  the  base 
pressure  becomes  infinite  as  the  normalized  height  V>  (hereafter  called 
simply  the  height)  tends  to  zero  which  results  in  an  infinitely  stiff 
spring  as  height  tends  to  zero.  Actually  the  base  pressure  cannot  exceed 


the  total  pressure,  Rj-  .  We  see  from  Figure  1  that  the  two  thick  jet 
theories  satisfy  this  condition.  The  base  pressure  also  approaches  the 
total  pressure  in  such  a  way  that  the  spring  stiffness  falls  to  zero  as 
the  height  approaches  zero.  The  spring  stiffness  for  small  height  de¬ 
pends  to  a  considerable  extent  on  the  derivative  of  the  normalized  momentum 
with  respect  to  the  normalized  base  pressure  Tb  ,  which,  as  shown 

d.  Po 

in  Figure  1,  tends  to  become  very  large  and  differs  considerably  for 
the  two  thick  jet  theories .  For  this  reason,  even  the  result  for  thick 
jet  theory  can  not  be  relied  upon  for  very  small  values  of  the  height 
(e.g.  V)  less  than  1.5 ).  We  see  from  these  Figures  that  the  lift  effi¬ 
ciency  factor,  VI  i_  ,  has  very  little  effect  on  the  value  of  the  natural 
frequency,  C*->^  ,  and  that  this  frequency  is  decreased  as  the  nozzle  to 

base  area  ratio,  CL  ,  is  increased.  Comparing  Figure  ^  with  Figure  5> 
shows  that  the  heave  frequency,  COwA  >  is  decreased  for  small  height 
and  increased  for  large  height  when  the  jet  inclination  angle,  P>0  ,  is 
changed  from  0  to  ^5  degrees .  The  nozzle  to  base  area  ratio,  Oo  ,  cannot 
vary  over  too  great  a  range  due  to  consideration  of  lift  efficiency.  (i.e., 
for  efficiency  the  vehicle  should  operate  at  small  values  of  height  to  dia¬ 
meter  ratio,  , where  D  is  an  equivalent  diameter  for  non-circular  GEMs, 

and  when  this  ratio  is  small  the  optimum  nozzle  to  base  area  ratio  is  small.) 

The  values  obtained  experimentally  (experimental  procedure  described  in 
Appendix  B)  are  shown  in  Figure  6.  Also  included  in  this  Figure  are  the 
appropriate  theoretical  curves  (  Qj=  ^  ,  fio  -  ^-S0  >  ~fc  =■  I  ) 

and  the  thin  jet  theory  expression  ^  is  the  acceleration  of 

gravity) . 
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The  theoretical  expression  for  the  damping  ratio  (43)  has  been  plotted  in 
Figure  J.  It  has  been  normalized  in  the  form  ^  (where  R_  is  the 

specific  lift  or  weight  per  unit  base  area  and  t  is  the  jet  thickness).  The 
nozzle  width,  t  ,  represents  the  effect  of  scaling  since 
a  function  of  the  acceleration  of  gravity,  the  air  density,  p  ,  and  non- 
dimensional  parameters.  Thus,  increasing  the  specific  lift  (i.e.,  ^  )  de¬ 
creases  the  damping  ratio,  ,  but  has  no  effect  if  the  size  of  the  machine 

is  increased  proportionally.  We  note  from  this  Figure  that  the  effect  of  the 
jet  to  base  area  ratio,  CL  ,  is  considerably  greater  than  that  of,  |3o  >  and, 

.  This  dependence  on  CL  is  physically  reasonable.  The  experimental  values 
obtained  on  the  eight  foot  model  are  shown  for  different  values  of  specific  lift, 

Pi_  ,  in  Figures  o,  9  and  10.  The  comparison  with  the  theoretical  curves  is 
fairly  good.  The  magnitude  is  appreciably  correct  and  there  is  no  noticeable 
dependence  on  height,  \o  .  Figure  11  shows  how  the  damping  ratio,  ,  varies, 

for  constant  RPM,  with  specific  lift,  PL  ,  as  a  parameter  along  the  curve. 

Finally,  in  Figure  12  we  have  plotted  a  number  of  calculated  and  experi¬ 
mental  values  of  the  roots  of  the  heave  stability  equation.  The  variation  with 
the  height,  Vl  ,  is  shown  for  several  values  of  specific  lift,  Pi_  . 

3.2.  Pitching  motion. 

In  an  earlier  section  of  this  report  momentum  considerations  were  utilized 
to  obtain  an  estimate  of  the  moments  acting  on  a  GEM.  In  particular,  an  ex¬ 
pression  was  obtained  for  the  slope  of  the  moment  versus  attitude  curve  for  a 
machine  trimmed  parallel  to  the  ground.  Referring  to  this  expression  (Equation  60), 
we  note  that  the  two  terms  in  the  numerator  depend  on  ('-«)  and  (\-fii.)  respectively, 
so  that  if  both  the  cross  flow  dissipation  factor,  C*  ,  and  the  lift  efficiency 
factor,  ,  are  equal  to  one,  the  moment  derivative  is  identically  zero.  The 

behavior  of  the  moment  curve  slope  as  a  function  of  these  parameters  will  be  com¬ 
pared  with  the  experimental  results  obtained  on  the  eight  foot  model. 


t. 


only 
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Recall  that  Equation  60  for  the  moment  coefficient  derivative. 


dOn 

ae  ' 

shows  that  when  &  =  I  (i.e.,  no  dissipation  of  cross  flow)  this  derivative 
is  always  negative.  Thus,  the  momentum  theory  predicts  that  in  the  absence 
of  jet  dissipation,  the  GEM  would  be  stable  at  all  altitudes.  We  also  see 
from  the  same  expression  that  if  the  lift  efficiency  factor,  i_  ,  were 
equal  to  one  (i.e.,  if  the  lift  were  actually  as  predicted  by  the  momentum 
theory)  then  the  GEM  would  be  unstable  at  all  altitudes.  Experience  has 
shown  that  typical  values  of  the  lift  efficiency  factor,  If) j_  ,  are  in  the 
range  0.7  to  0.8.  (Reference  27,  for  example,  suggests  a  value  0.8  for  well 
designed  machines.)  Since  all  machines  tested  to  date  are  stable  at  a  very 
low  altitude  and  become  unstable  as  the  altitude  is  increased,  the  cross  flow 
dissipation  factor,  cx  ,  must  be  less  than  one  for  the  moment  slope  to  give 
results  which  correspond  with  experiment .  Figures  13  and  11+  show  how  curves 
of  the  moment  coefficient  derivative,  versus  height  are  affected  by 

changes  in  the  cross  flow  dissipation  factor,  (X  .  Larger  values  of  (X  are 
not  included  since  they  result  in  unreasonably  high  stability  crossover  points 
(i.e.,  values  of  height,  ,  where  the  vehicle  becomes  neutrally  stable).  We 
see  that  the  curves  differ  considerably  even  for  small  changes  in  CK  .  In¬ 
creasing  (X  ,  decreases  the  maximum  of  the  moment  coefficient  derivative  and 
increases  the  crossover  point .  Figure  14  which  is  plotted  for  a  larger  value 
of  the  lift  efficiency  factor,  Q.  u  ,  shows  the  same  general  behavior  with 
respect  to  the  cross  flow  dissipation  factor,  Q(  ,  although  the  crossover 
point  does  not  increase  as  rapidly.  Figure  17  in  which  the  jet  inclination 
angle,  $0  ,  is  zero,  also  shows  the  same  behavior  as  a  function  of  the 

cross  flow  dissipation  factor. 
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The  behavior  of  the  curves  for  different  values  of  the  lift  efficiency 
factor,  I7l.  ,  is  shown  in  Figures  18  and  19.  These  two  sets  of  curves 
are  drawn  only  for  the  cross  flow  dissipation  factor  0(  «*  O  .  The  general 
character  of  the  curves,  however,  is  the  same  for  all  values  of  this  factor. 
We  note  that  the  behavior  of  the  moment  coefficient  derivative, 
with  varying  height  is  quite  different  for  the  different  values  of  the 
lift  efficiency  factor,  .  As  was  pointed  out  previously,  repre¬ 

sents  the  fact  that  the  actual  lift  produced  by  the  base  pressure  falls  short 
of  that  computed  on  the  basis  of  a  uniform  base  pressure.  There  are  probably 
a  number  of  phenomena  not  included  in  our  simplified  approach  which  could 
account  for  this  decrease  in  lift.  Probably  the  most  significant  of  these 
is  the  vortices.  Smoke  studies  of  the  flow  under  the  base  show  that  there 
are  vortices  near  the  jets  under  the  base  which  reduce  the  pressure  on  the 
base  below  the  pressure  which  acts  near  the  center  of  the  base. 

When  the  lift  efficiency  factor,  He  ?  '  the  GEM  is  unstable  for  all 
values  of  height,  V)  .  For  smaller  values  of  the  lift  efficiency  factor, 

Hl  ,  the  machine  becomes  stable  for  some  heights,  but  the  cross  over 


cL  O-i 
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height  and  the  maximum  unstable  value  of  the  moment  coefficient  derivative. 


varies  considerably  for  different  values  of  the  lift  efficiency 

cL  © 

factor,  L  .On  comparing  Figures  18  and  19  we  see  that,  for  a  particu¬ 


lar  value  of  the  left  efficiency  factor,  the  machine  becomes  stable  only 
at  a  much  lower  height,  Vi  ,  when  the  jet  inclination  angle,  (30,  equals 
zero  and  that  the  maximum  unstable  moment  slope  for  this  case  is  also  con¬ 


siderably  greater. 


53 


varies 


Figure  1 6  shows  how  the  moment  coefficient  derivative 


dCei 


for  fixed  values  of  lift  efficiency  He.  and  jet  to  base  area  ratio,  CL  , 
as  the  jet  inclination  angle,  (3  o  ,  varies  from  zero  degrees  to  60  degrees 
We  see  that  the  GEM  becomes  increasingly  stable  as  this  angle  increases. 
From  Figure  15  we  see  that  the  jet  to  base  area  ratio,  CL  f  also  has  some 

influence  on  the  value  of  the  moment  coefficient  derivative  . 

a-  © 

This  influence,  however,  is  much  less  than  the  influence  of  other  parameters 
We  now  consider  the  experimental  results.  Static  moments  were  measured 
for  different  values  of  pitch  (or  roll)  angle,  ©  ,  at  a  number  of  height 
and  power  settings  throughout  the  range  available  in  the  experimental  set  up 
Some  typical  curves  obtained  in  this  manner  are  shown  in  Figure  21.  For  the 
lower  heights  the  moment  curve  is  very  nearly  a  straight  line  and  has  a 
stable  slope.  As  the  height  is  increased,  the  moments  reach  a  maximum  at 
some  angle  and  then  drop  off.  The  angle  at  which  this  maximum  occurs  be- 
comes  larger  as  the  height  is  increased.  The  moment  curves  were,  however, 
nearly  linear  near  zero  angle  for  all  heights. 

It  will  be  noticed  that  in  most  of  the  curves  zero  moment  does  not 
occur  exactly  at  zero  angle.  This  can  be  assumed  to  be  due  to  slight 


asymmetries  in  the  flow  field  of  the  model.  The  slopes  of  the  moment 
curves  at  zero  incidence  are  plotted  as  a  function  of  height,  V~>  ,  in 

Figure  20  along  with  the  appropriate  theoretical  curves.  A  number  of 
curves  were  computed  to  determine  which  combination  of  the  cross  flow 
dissipation  factor,  0(  ,  and  the  lift  efficiency  factor,  Qi_,  would 
give  the  best  match  to  the  experimental  data.  It  was  found  that  the 
closest  matches  were  obtained  for  H.  i_  between  0.7  and  0.8  and  for  C>( 
between  0  and  0.2.  All  of  the  theoretical  curves  have  a  maximum  for 


a  height,  Vi  ,  of  the  order  of  10,  whereas  the  experimental  points  seem  to 

be  increasing  with  height  even  at  equal  to  11.  We  should  note  here, 

however,  that  one  of  the  conditions  for  the  theory  to  be  valid  is  that  the 

L_ 

height  to  diameter  ratio,  ,  not  be  too  large,  since  for  sufficiently 

L, 

large  -Ih_  the  jets  on  the  right  and  left  of  the  two  dimensional  model 
may  come  together.  In  this  case,  the  moment  behavior  may  be  quite  dif¬ 
ferent  from  that  shown  by  the  present  theory. 

For  small  values  of  the  height,  W  ,  the  theoretical  curves  drop  off 
much  more  rapidly  than  the  experimental  curve .  This  is  probably  due  to  the 
fact  that  one  of  the  assumptions  of  the  analysis  becomes  weak  when  the 
height  becomes  sufficiently  small.  This  is  the  assumption  that  the  total 
pressure  in  the  chamber,  Pj-  ,  could  be  considered  constant.  Also,  the 
derivative  of  the  normalized  jet  momentum  with  respect  to  the  normalized 
base  pressure,  ^~ta  ,  as  shown  in  Figure  1,  becomes  very  large  at  very 

d-PB 


small  height,  W  ,  which  is  probably  unrealistic.  These  results  seem  to 
indicate  that  the  theory  is  satisfactory  for  estimating  the  magnitude  of 
the  slope  of  the  moment  curve,  as  veil  as  its  behavior  with  height  changes 
and  changes  in  the  various  physical  parameters  of  the  vehicle.  Accurate 
values  of  the  moment  slope  can  not  be  determined,  however,  due  to  the  sensi¬ 
tivity  of  the  results  to  the  lift  efficiency,  ,  and  the  dissipation 

factor  0(  •  Although  these  parameters  can  be  estimated,  they  can  not  be 

easily  determined  accurately.  Also,  to  improve  the  accuracy  of  the  results 
for  small  height,  w  ,  it  would  probably  be  necessary  to  refine  the  theory 
and  include  the  effects  of  changes  in  total  pressure  and  vortices. 
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Using  the  present  theoretical  approach  it  was  not  possible  to  find  a 
satisfactory  explanation  for  damping  in  the  pitching  motion.  Experiments 
seemed  to  indicate,  however,  that  there  was  a  considerable  amount  of  damp¬ 
ing  in  this  mode  of  motion.  Results  of  these  experiments  are  shown  in 
Figure  22.  These  curves  indicate  that  the  damping  decreases  quite  rapidly 
as  height,  Vi  ,  is  increased.  The  damping  itself,  however  was  relatively 
high  in  all  of  the  tests . 

These  results  indicate  that,  for  all  but  the  very  lowest  heights 
tested,  the  characteristic  roots  of  the  pitching  mode  are  both  real,  one 
lying  to  the  left  of  the  imaginary  axis  and  one  to  the  right.  If  there  had 
been  no  damping  in  the  system  the  two  roots  would  have  been  placed  symmetri¬ 
cally  with  respect  to  the  imaginary  axis.  The  damping  moves  both  roots  to 
the  left.  In  Figure  23  we  see  how  these  roots  move  along  the  real  axis  as 
the  altitude  is  varied.  The  abscissa  of  the  curve  shown  represents  the 
position  of  the  two  roots  along  the  real  axis  for  a  given  height.  As  the 
height,  V)  ,  decreases,  the  divergent  root  moves  toward  the  origin,  pass¬ 
ing  through  the  value  zero  at  =  3-3>  while  the  convergent  root  at  first 
moves  away  and  then  starts  back  toward  the  origin  at  about  the  same  height, 
r)  at  which  the  divergent  root  passes  through  the  value  zero.  The  two 
roots  come  together  at  a  position  -7.3  on  the  real  axis  for  \r\  equal  to  1.4. 
For  heights  lower  than  this  the  roots  would  represent  a  damped  oscillation. 
This  last  point  is  the  result  of  extrapolation. 
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4.  FORWARD  FLIGHT  CONSIDERATIONS 


In  considering  the  hovering  GEM  we  were  able  to  present  a  theory  which 
predicted  the  most  important  dynamic  characteristics  of  the  machine's  behavior. 
We  were  also  able  to  present  experimental  results  with  which  to  check  these 
predictions.  For  the  GEM  in  forward  flight  the  situation  is  considerably  more 
difficult.  No  adequate  theory  is  available  and  the  small  amount  of  experi¬ 
mental  results  available  are  for  fairly  high  free  stream  velocities  (greater 
than  30  feet  per  second).  For  this  reason,  the  discussion  given  here  will  be 
confined  to  summarizing  what  is  known  about  the  forward  flight  regime  and  to 
making  a  few  observations  as  to  how  this  might  effect  the  dynamic  character¬ 
istics  of  the  GEM  in  low  speed  flight. 

The  character  of  the  flow  around  the  GEM  changes  considerably  as  the  free 
stream  velocity  is  increased.  Depending  on  the  dynamic  pressure,  the  flow 
field  will  be  similar  to  one  of  forms  illustrated  as  A  and  B  below. 


A  Low  Velocity  B  High  Velocity 


If  the  freestream  velocity  is  zero  there  is  no  cross  flow  when  the  GEM  is  parallel 
to  the  ground.  When  the  free  stream  velocity  is  not  zero,  the  static  pressure  in 
front  of  the  jet  will  be  increased.  The  difference  between  the  base  pressure  and 
the  ambient  pressure  becomes  less  so  that,  if  the  height  remains  constant,  the 
forward  jet  need  not  be  as  strong.  Thus,  part  of  the  jet  will  pass  under  the 
GEM  as  shown  in  A  above.  As  the  dynamic  pressure  increases,  the  strength  of  the 


# 


cross  flow  will  increase  until  finally  all  of.  the  jet  passes  under  the  base  as 
shown  in  B  above.  The  velocity  at  which  the  flow  pattern  A  changes  to  pattern 
B  will  be  approximately  determined  by  equating  the  free  stream  dynamic  pressure 
to  the  base  pressure. 
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for  small  l-»  .  This  will  be  only  a  very  rough  estimate  since  the  vortex 

ahead  of  the  jet  (shown  in  A)  will  influence  the  pressure  there.  We  neglect 
this  effect.  Observe  that  the  velocity  at  which  the  flow  change  occurs  will 
increase  with  increasing  specific  load.  Smoke  studies  described  in  Reference  29 
seem  to  indicate  that  the  jet  turns  back  under  the  base  at  a  dynamic  pressure 
somewhat  lower  thur.  the  one  corresponding  to  the  base  pressure .  If,  however, 
the  specific  load,  ,  is  very  low,  flow  pattern  B  might  be  reached  before 

the  dynamic  pressure  is  high  enough  to  make  control  surfaces  effective .  We 
shall  consider  here  only  forward  flight  at  low  speeds  where  control  surfaces 
are  not  effective.  The  speed  at  which  the  control  surface  becomes  useable 
would  depend  on  the  particular  design  under  consideration. 

In  hovering  we  could,  with  sufficient  accuracy,  consider  the  motion  to 
have  only  two  degrees  of  freedom  (i.e.,  forces  introduced  by  horizontal  motion 
about  the  zero  velocity  condition  could  be  neglected),  we  now  have  three  de¬ 
grees  of  freedom.  For  low  speed  flight  of  a  level  GEM,  however,  the  heave 
degree  of  freedom  is  nearly  uncoupled  from  the  pitch  and  forward  velocity  de¬ 
grees  of  freedom.  There  is  a  small  coupling  because  there  is  some  lift  change 
with  either  a  change  in  forward  velocity  or  a  change  in  attitude.  The  former 
lift  change  should  be  small  if  the  forward  velocity  is  sufficiently  small,  and 
the  latter  lift  change  should  be  zero  (or  nearly  so)  if  the  aircraft  is  initially 
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trimmed  in  a  level  attitude.  Thus,  in  order  to  keep  the  analysis  simple,  we 
neglect  the  heave  degree  of  freedom  in  forward  flight  and  consider  only  the 
attitude  and  translational  degrees  of  freedom.  This  does  not  imply  that  the 
variation  of  vertical  force  with  changes  in  the  steady  state  velocity  will  be 
small.  However,  such  variation  should  introduce  no  severe  dynamic  problems. 

The  perturbation  equations  describing  the  GEM  in  forward  flight  will 
thus  be, 
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The  rate  term,  0  ,  which  arises  in  conventional  aircraft  primarily 

because  of  a  change  in  the  angle  of  attack  of  the  horizontal  tail,  can  be 
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neglected  at  low  velocities.  We  then  need  to  estimate 


moment  change  with  forward  velocity; 


the  pitching 


a© 


,  the  horizontal  force  change  with 


pitch  angle;  ,  the  horizontal  force  change  with  velocity;  ,  the 

ov  <3© 

moment  due  to  attitude  change;  and  ,  the  damping  moment.  We  would  like 

to  know  how  these  derivatives  change  with  changes  in  the  height,  Vi  ,  and  the 
forward  velocity,  V 

The  forces  acting  on  the  GEM  in  the  horizontal  direction  are: 

a.  the  thrust  from  auxiliary  propulsion 

b.  the  component  of  lift  vector  in  horizontal  direction 

c .  the  momentum  drag 

d.  the  parasite  drag  (skin  friction  and  form  drag) 

e.  the  interference  thrust  (momentum  recovery). 
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Additional  drag  will  be  encountered  in  flight  over  water;  however,  this 
effect  will  not  be  considered  here.  For  a  more  complete  discussion  of  drag 
forces  acting  on  GEMs,  see  Reference  21  (Andu  and  Miyashita). 

Thus,  (with  the  X-axis  parallel  to  the  ground)  the  horizontal  force  can 
be  expressed: 

'X  =  ~T  Ctyl  ©  —  ©  -yu  V  -  V2  Cop  +  TVn£  (78) 

where  T  =  propulsion  thrust 

L  =  lift 

JW  =  mass  flow 

C*C>p  =  parasite  drag  coefficient 
T>nX  =  interference  thrust 
We  shall  assume  that  the  thrust  is  constant  and  that 
O  =  0  in  equilibrium;  thus  “towr*  ©  =  © 

W  (machine  weight) 

Cop  is  approximately  independent  of  both  V  and  © 

Thus  we  obtain. 
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Little  is  known  about  the  interference  thrust  .  The  following 

experimental  curves  from  Reference  26  (Sweeney  and  Nixon),  however,  show  that 


it  is  strongly  dependent  on  velocity. 
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Andu  and  Miyashita  (Reference  21 )  have  suggested  the  following  mechanisms 
for  this  term, 

a.  the  disappearance  of  dead  air  at  the  air  intake  and  nozzle  exit 
(which  was  present  in  power-off  parasite  drag  experiments),  and  the  re¬ 
duction  of  friction  drag  on  the  base,  reduce  the  parasite  drag  in  a  power- 
on  condition; 

b .  the  ingested  and  ejected  air  bring  about  a  favorable  effect  on  the 
external  flow; 

c.  the  ejected  air  has  energy  in  the  form  of  total  jjressure,  some  part 
of  which  is  converted  to  propulsive  work  in  the  process  of  expansion. 

Because  of  the  present  lack  of  understanding,  the  derivatives  of  this  term 
will  have  to  be  estimated  from  experimental  data.  From  the  curve  shown  above, 
it  appears  that  the  interference  thrust  decreases  with  increasing  velocity.  No 
data  is,  however,  available  to  estimate  its  behavior  with  pitch  angle. 
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There  will  be  moments  acting  on  the  machine  for  the  following  reasons: 

a.  change  in  the  base  pressure  distribution 

b.  change  in  pressure  distribution  over  the  exterior  of  the  machine 

c .  change  in  jet  strength  at  the  front  and  back  of  the  machine . 

As  soon  as  the  GEM  starts  to  move  forward,  the  pressure  distribution  over  the 
upper  surface  and  around  the  exterior  of  the  jet  will  change .  One  would  at 
first  think  that  the  increased  pressure  at  the  front  of  the  machine  would 
cause  the  base  pressure  just  inside  of  the  front  jet  to  increase  and  hence 
give  a  positive  nose-up  pitching  moment .  Experiments,  however,  have  shown 
that,  for  low  velocities  (corresponding  to  condition  A  above),  this  is  not 
the  case  (see  Reference  19).  This  effect  would  appear  to  be  attributable  to 
the  presence  of  the  large  vortex  at  the  front  of  the  machine .  Measurements 
of  the  base  pressure  distribution  given  in  Reference  19  show  that  the  pressure 
is  lowest  at  the  front  of  the  machine  and  increases  monotonically  (although  not 
always  linearly)  to  the  rear  of  the  machine.  This  results  in  a  negative  pitch¬ 
ing  moment .  Because  of  the  presence  of  the  forward  vortex  and  the  fact  that 
air  is  being  drawn  into  the  top  of  the  machine,  it  would  be  very  difficult  to 
make  any  general  statements  about  the  pressure  distribution  over  the  upper  sur¬ 
face  of  the  machine.  It  will  certainly  be  considerably  different  from  the 
pressure  distribution  of  an  unpowered  machine  and  it  will  probably  depend 
strongly  upon  the  particular  exterior  shape  of  the  GEM. 

Moments  arising  directly  from  the  change  in  jet  strength  will  be  small  as 
long  as  the  jet  area  is  small  compared  to  the  base  area.  If  we  neglect  this 
contribution  of  the  jets  we  can  then  express  the  moment  as  follows: 

M=  Mext  -1-  -3-r2  aP  (80) 


62 


where 
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M  ext  =  moment  due  to  the  pressure  changes  on  the  upper  surface  of  the  machine 

AP  =  difference  between  the  base  pressure  at  the  front  and  the  back  of 

the  machine 

=  equivalent  radius  of  the  machine 

=  a  factor  which  accounts  for  the  three  dimensional  character  of  the 
machine,  as  well  as  for  differences  between  the  actual  and  a 
linear  pressure  distribution. 

We  are  mainly  interested  in  estimating  the  moment  change  with  change  of  atti- 
tude  angle,  ?.r\  >  and  the  moment  change  due  to  change  in  forward  velocity  . 

0c7  Q V 

We  can  express  these  as  follows : 
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The  damping  in  pitch,  «  ,  will  probably  not  be  much  effected  by  the 

small  forward  velocities  and  will  be  considered  to  be  approximately  the  same 
as  for  a  hovering  GEM. 

Because  of  the  complexity  of  the  flow  under  the  base  and  across  the  surface 
of  the  GEM,  the  two  derivatives  and  must  be  determined  for  any  specific 

d  ©  ft V 

machine  by  experimental  means.  Using  the  results  of  some  wind  tunnel  tests  which 
have  been  made,  however,  we  can  make  some  general  qualitative  observations  con¬ 
cerning  the  manner  in  which  these  two  derivatives  can  be  expected  to  behave. 

Consider  first  the  moment  change  due  to  velocity,  •  We  can  obtain 

some  feeling  for  the  behavior  of  this  derivative  from  the  results  of  Reference  19 • 
The  curves  below  show  the  qualitative  behavior  of  the  moment  coefficient,  0 , 
with  changes  in  velocity  for  different  heights.  (The  notation  is  ours  and  not 
that  of  the  Reference  cited). 
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The  parameter  0/J  is  defined  to  be: 


c  *  TJl 
^  Vaf>s/2  5 

where  J  is  jet  strength  per  unit  le'ngth;  is  the  circumference;  and  S 

is  the  base  area  of  the  GEM. 

If  we  equate  the  dynamic  pressure  to  the  base  pressure  and  use  the  above 
definition  of  we  find  that. 


l  rU  l 


4  V 


In  the  curves  shown  above,  the  sharp  changes  occur  at  a  value  of  -p=r — 

(-/J- 

approximately  equal  to  the  one  given  in  this  Equation.  Thus,  the  region  of 
interest  to  us  in  these  curves  is  only  that  small  portion  lying  to  the  left  of 
the  sharp  slope  changes  (i.e.,  to  the  left  of  the  dotted  line  in  the  Figure). 


We  note  from  the  curves  that,  for  very  small  height  to  diameter  ratio,  Vb  > 

...  •  ,  ,  j_ i_ l  „ _j_  i -|  - 


the  derivative 


£>V 


is  positive  and  that,  for  the  larger  values  of 


it  is 


negative  for  small  forward  velocities  and  changes  sign  rapidly  at  some  velocities, 
As  Wd  increases  the  magnitude  of  the  slope,  ,  becomes  greater.  This 

would  seem  to  indicate  that  an  accelerating  GEM  would  experience  rapid  changes 
in  dynamic  characteristics  if  it  were  operating  at  sufficiently  large  value 


of 


Vd  • 


We  now  consider  the  derivative  A-  •  Almost  no  information  is  available 

2>© 

concerning  the  effect  of  forward  velocity  on  this  derivative  in  the  region  in 
which  we  are  interested  (i.e.,  bp'/2  <  Pa  )•  This  is  because  of  the  fact 
that  all  moment  information  which  is  available  comes  from  wind  tunnel  tests 
which  cannot  be  conducted  with  accuracy  at  less  than  about  30  feet  per  second. 

In  order  to  conduct  experiments  at  lower  velocities  a  facility  such  as  the 
Princeton  University  Forward  Flight  Facility  would  be  required.  Although  some 
ground  effect  experiments  have  been  conducted  in  this  Facility  with  a  four  foot 
model,  none  have  yet  included  pitching  moment  measurements. 

Sweeney  and  Nixon  (Reference  30)  have  stated  that  the  Princeton  20-foot  GEM 
(P-GEM)  becomes  statically  stable  in  forward  flight,  although  it  is  quite  un¬ 
stable  in  hover  (at  full  power  with  stabilizing  slots  closed).  They  remark, 
however,  that  this  might  be  associated  with  the  fact  that  the  P-GEM  achieves 
forward  propulsion  by  tilting  the  nose  down,  and  that  the  stability  increase 
could  be  associated  either  with  the  slight  height  loss  in  this  attitude  or  to 
the  attitude  itself.  Further  tests  with  a  leveled  machine  will  probably  be  re¬ 
quired  to  determine  whether  the  note  of  optimism  in  these  results  is  justified. 
However,  one  might  conclude  from  these  results  that  the  effect  of  low  forward 
velocities  will  at  least  not  decrease  the  static  stability. 
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5.  PRELIMINARY  CONTROL  SYNTHESIS* 

5.1.  Stability  augmentation 

General  experience  with  ground  effect  machines,  as  well  as  the  theoretical 
results  presented  in  the  first  part  of  this  report,  indicate  consistently  that 
a  peripheral  jet  machine  exhibits  static  instability  in  attitude  well  within 
the  height  range  where  it's  performance  is  still  economical.  Stabilization 
by  aerodynamic  means,  such  as  slots,  has  to  be  paid  for  in  performance. 
Artificial  stability  augmentation  by  means  of  automatic  feedback  control 
suggests  itself  as  a  powerful  and  convenient  alternative. 

The  development  of  the  ideas  concerning  artificial  s  .ability  augmentation 
consists  of  three  steps: 

a.  description  of  the  uncontrolled  system 

b.  establishing  design  principles  for  the  controlled  system 

c.  conclusions  to  be  drawn  from  (a)  and  (b)  as  to  the  design  of  a 
controller,  as  well  as  consideration  of  the  limitations  of  its  use. 

In  the  following  we  shall  consider  first  pitch  stabilization  in  hover,  and 
then  discuss  the  effect  of  forward  velocity  on  this  proposed  stabilisation. 

5.1.1.  Description  of  the  uncontrolled  system.  In  the  initial  approach  we 
are  going  to  consider  only  the  attitude  degree  of  freedom  (pitch  or  roll).  In 
the  previous  discussion  of  the  equations  of  motion,  the  attitude  degree  of 
freedom  has  been  singled  out  as  causing  the  most  significant  stability  problems. 
Tests  and  theory  have  shown  that  the  aerodynamic  "spring"  of  this  second  order 
system  varies  from  positive  (stabilizing)  to  negative  (destabilizing)  values. 
With  increased  hovering  height  the  negative  "spring  constant"  reaches  a  rather 
flat  maximum  absolute  value.  The  practically  useful  performance  range  of  the 
GEM  reaches  into  the  region  of  this  maximum. 

*Methods  used  below  are  discussed  in  servomechanism  texts  (e.g.  Reference  31) 
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The  transfer  function  for  one  degree  of  freedom  in  attitude  can  be  written 


in  the  following  form: 


G 


e 


Ac/t 

Cs-s,)(.s-  Si) 


where 

s,  -  coP  (  -  r+  +  ?*) 


^2  -  ^ p  ("  ) 


I 


(83) 


With  no  damping,  ^=o  >  S(=  - Sa  =  COp  .  The  effect  of  the 

damping  on  the  location  of  the  roots  of  a  second  order  system  with  a  negative 
spring  is  shown  in  Figure  25.  Our  test  results  (see  Figures  22  and  23)  have 
shown  the  existence  of  fairly  high  damping  throughout  the  useful  height  range, 
although  no  theoretical  explanation  for  this  phenomenon  could  be  found. 

It  seems  reasonable  to  assume  that  the  general  behavior  of  the  roots 
(  S,  and  S2.  )  would  be  similar  to  the  behavior  indicated  in  Figure  23,  and 
that  this  would  be  independent  of  the  size  of  the  machine. 

The  effect  of  scaling  on  the  "undamped  natural  frequency",  COp  ,  is  the 
following : 


5M  /  ae 
I 


M  ^ 


(84) 


where  is  the  radius  of  gyration  and  -=-  is  the  form  in  which  the 

moment  derivatives  are  plotted  (Figures  13  through  20).  With  and 

f 

\N=  L  we  ob  tain 
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7.  Or  /  r  \a  (  '  5Cm  \ 

=  ^  -b  (  p^./  \  r  50  / 

This  expression  indicates  a  decrease  in  the  natural  frequency,  COp  , 
inversely  proportional  to  the  square  root  of  the  size  of  the  machine,  repre¬ 
sented  in  this  equation  by  the  nozzle  width,  "t  .  It  should  be  noted  that  the 
explicit  dimension,  "t  ,  in  this  expression  reflects  the  size  of  the  machine; 
the  effect  of  the  nozzle  width  on  the  moment  derivative  has  been  indicated  in 
Figure  15  in  terms  of  the  nozzle  area  to  base  area  ratio. 

In  the  damped  case  a  change  in  the  natural  frequency,  COp  ,  means  a  pro¬ 
portional  change  in  the  distance,  (Op  +  ,  of  the  roots  from  their 

arithmetic  mean,  —  COp  .  The  effect  of  a  change  in  the  damping  ratio  on 
the  arithmetic  mean  is  also  proportional.  Its  effect  on  the  distance  between 
the  roots  is,  however,  less  than  proportional. 

Because  of  the  simplifying  assumptions  made  in  the  derivation  of  the 
moment  derivative,  and  because  of  the  lack  of  knowledge  about  the  damping, 
we  cannot  expect  highly  accurate^  theoretical  results.  The  use  of  the  theo¬ 
retical  results  together  with  the  experimental  results,  however,  yields  suffi¬ 
cient  information  for  an  approximate  model  to  be  used  in  a  preliminary  control 
system  synthesis.  We  observe  (see  Figure  23)  that  the  root  on  the  left  hand 
side  does  not  change  much  throughout  the  useful  range  of  heights.  For  practical 
purposes  this  root  can  be  assumed  to  be  constant.  It  is  suggested  that  its 
value  could  be  determined  by  assuming  a  certain  damping  ratio,  ,  at  the 

height  where  is  maximum:  =  ^Pmax  V*  )  •  The  other 

root  would  be  determined  from  the  fact  that  the  product  of  the  roots  is 

z 

S,  Si  —  Co p  .  Either  theoretical  or  experimental  moment  derivatives 
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could  be  used  to  determine  the  natural  frequency,  COp  .  In  order  to  approxi¬ 
mate  the  mean  value  of  the  root,  St  ,  a  correction  factor  should  be  used, 

Sim  —  '•  2.  Si  in  .  Assuming  ^  ^  =  1  and  using  the  theoretical  curve 

(Figure  20)  for  the  values  =  O  we  obtain:  GOp (rva^  =  -4 .  (o 

at  =  ST  and  S^—  43-3  .  This  value  and  the  curve  calculated  for  S2  are 

shown  as  dashed  lines  in  Figure  23.  The  effect  of  deviations  from  our  somewhat 
arbitrary  approximations  will  be  discussed  later. 

We  return  now  to  the  expression  (Equation  83)  for  the  transfer  function 
describing  the  uncoupled  attitude  response  to  a  control  input,  Se  •  The  gain, 

Ac  >  is  actually  an  approximation  to  the  transfer  function  which  would  accu¬ 
rately  express  the  relationship  between  the  control  moment  and  the  control  input. 
This  relationship  varies  with  the  height  and  generally  may  involve  some  dynamics. 
For  example,  if  throttling  of  a  portion  of  the  jet,  or  inclination  of  the  jet,  is 
used  in  order  to  produce  control  moments,  a  part  of  the  moment  develops  immedi¬ 
ately  in  proportion  to  the  change  in  jet  momentum.  Another  part  of  the  control 
moment,  however,  will  develop  only  after  the  flow  pattern  under  the  base  has 
changed.  At  the  present  time,  no  information  is  available  on  the  magnitude  and 
the  lag  of  this  part  of  the  control  moment  which  must  depend  largely  upon  the 
height  above  the  ground  as  well  as  size  of  the  machine.  This  is  one  reason  why 
a  considerable  margin  must  be  introduced  in  the  design  criteria. 

5.1.2.  Design  principles  for  the  controlled  system.  Establishing  reasonable 
and  realistic  performance  criteria  for  the  feedback  control  of  a  vehicle  so 
vaguely  understood  as  the  ground  effect  machine,  is  perhaps  the  most  intricate 
step  in  the  entire  synthesis  procedure.  The  following  discussion  is  actually  not 
restricted  to  ground  effect  machines  but  is  of  broader  validity. 

Our  design  criteria  must  be  based  on  criteria  for  controllability  by  a  pilot. 
In  recent  years  considerable  progress  has  been  made  toward  a  technically  meaningful 
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formulation  of  the  requirements  which  have  to  be  met  in  order  to  permit  a 
human  pilot  to  perform  satisfactory  control.  (See  References  22,  23  and  24.) 

This  approach  is  based  on  the  fact  that  the  human  pilot  is  actually  an  element 
in  a  feedback  loop  closed  around  the  controlled  airframe  or  vehicle. 

Extensive  experiments  have  shown  that  the  response  of  a  human  operator 
in  a  feedback  loop  can  be  considered  to  be  made  up  of  two  parts.  One  part 
is  linearly  correlated  with  the  pilot's  input;  the  remainder,  usually  called 
the  remnant,  cannot  be  linearly  correlated  with  the  input.  The  latter  part 
is  presumably  attributable  to  nonlinearities  in  the  pilot's  action,  to  noise 
introduced  by  the  pilot,  and  to  the  pilot's  response  to  inputs  other  than  the 
ones  considered  in  the  experiment.  The  following  basic  ideas  have  emerged  for 
describing  the  human  operator. 

a.  In  order  to  obtain  the  linearly  correlated  part  of  the  pilot's  res¬ 
ponse,  the  pilot  can  be  replaced  by  a  transfer  function  of  the-  form  (see 

/\p  <2  I  Ti_S 

I  +  S  \  +  Tx  S 

-reaction  time  delay  >  basically  not  adjustable 
=  neuromuscular  lag  by  the  pilot 

=  the  pilot's  equalization  characteristic 

=  the  pilot ' s  gain 

This  transfer  function  is  a  good  approximation  for  a  wide  variety  of 
control  tasks . 

b.  The  human  pilot  is  highly  adaptive.  He  is  able  to  change  his  gain,  AP, 
his  lead  time  constant,  Ti_  ,  and  his  lag  time  constant,  Tx  ,  over  a  wide  range. 


Reference  24): 
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numerous  experiments  indicate  that  his  efforts  are  aimed  at  providing  a  suf¬ 
ficient  phase  margin  (70-80  degrees )  for  the  open  loop  transfer  function  so  as 
to  assure  satisfactory  closed  loop  control. 

c.  The  pilot's  opinion  of  the  handling  qualities  of  the  machine  is  strongly 
correlated  with  the  gain  as  well  as  the  lead  and  lag  time  constants  he  must  use 
in  order  to  achieve  a  satisfactory  degree  of  stability.  This  opinion  deteriorates 
rapidly  when  he  reaches  the  limits  of  his  capabilities . 

Some  important  deductions  from  the  correlation  between  the  pilot's 
opinion  and  the  transfer  function  of  the  human  operator  are  the  following: 

a.  The  human  pilot  seems  to  have  the  best  opinion  of  tasks  in  which  only 
a  pure  gain  and  no  lead  or  lag  has  to  be  applied  on  his  part . 

b.  There  is  a  rather  flat  range  of  optimum  values  for  the  gain,  Ap  ,  of 
the  pilot.  A  considerable  range  (e.g.,  1:8  in  Reference  23)  is  acceptable. 

Too  low  or  too  high  a  gain  are  respectively  described  as  a  too  sluggish  res¬ 
ponse  or  a  too  sensitive  control. 

c.  The  pilot's  opinion  is  degraded  if  an  increasing  amount  of  equalization, 
especially  lead,  is  necessary.  Also  the  remnant  (i.e.,  uncorrelated  part  of  the 
pilot's  response)  is  relatively  small  if  little  equalization  is  needed,  and  it 
increases  with  an  increasing  effort  on  the  part  of  the  pilot  to  provide  "good" 
control . 

d.  The  bandwidth  of  the  closed  loop  determines  the  disturbance  (gust)  band¬ 
width  for  which  the  vehicle  can  be  kept  well  under  control.  Therefore,  the  pilot 
will  tend  to  adjust  his  gain  and  compensation  so  that  he  can  meet  this  requirement 
under  given  gust  conditions . 

These  extremely  simple  and  reasonable  principles  have  been  chosen  by  the 
authors  as  a  guide  in  determining  design  criteria  for  the  feedback  stabilization 


71 


of  the  GEM.  No  particular  quantitative  pilot  transfer  function  can  be  identi¬ 
fied  because  only  extensive  flight  experiments  with  ground  effect  machines  involv¬ 
ing  the  evaluation  of  pilots'  opinions  of  the  handling  qualities  could  justify 
such  a  cnoice.  We  can,  however,  utilize  the  principles  which  have  been  outlined 
without  being  specific  about  the  pilot's  transfer  function.  This  will  be  shown 
later . 

It  is  important  to  appreciate  that  the  approximations  involved  in  the  des¬ 
cription  of  the  dynamics  of  the  GEM  leave  a  considerable  amount  of  uncertainty 
in  the  synthesis  of  feedback  stabilization.  There  are  two  possible  choices  in 
dealing  with  this  problem.  The  difference  in  the  principle  of  these  two  approaches 
is  of  the  utmost  practical  importance.  The  nature  of  each  possible  choice  is  dis¬ 
cussed  below. 

a.  Adaptive  control:  We  may  specify  some  dynamic  performance  criterion  for 
the  vehicle  which  is  to  be  controlled  by  the  pilot.  This  leads  to  the  application 
of  the  adaptive  control  principle.  In  this  case,  we  make  up  our  mind  about  a 
feasible  optimum  vehicle  performance  and  keep  checking  the  responses  of  the  ve¬ 
hicle  to  its  inputs.  These  are  then  continuously  compared  with  the  performance 
criterion.  Adjustments  are  automatically  made  to  the  gains  of  an  inner  feedback 
loop  closed  around  the  vehicle,  as  well  as  possibly  in  an  equalizing  network  be¬ 
tween  the  pilot  and  the  inner  loop . 

Applying  the  adaptive  control  principle  relieves  us  of  the  burden  of  finding 
how  significant  stability  derivatives  vary.  The  price  we  have  to  pay  is  in  the 
increased  complexity  of  the  feedback  system.  Also,  we  must  face  some  possibly 
severe  problems  in  connection  with  the  dynamic  stability  of  a  fast  acting  adapt¬ 
ive  loop  closed  around  an  inherently  unstable  vehicle . 

b.  Nonlinear  compensation^  If  we  know  how  the  significant  stability  deri¬ 
vatives  vary,  we  can  make  the  adjustments  in  the  feedback  loop  on  the  basis  of 
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measuring  or  estimating  the  parameters  which  affect  the  derivatives.  If  the 
parameters  are  very  weakly  or  not  at  all  coupled  dynamically  to  the  controlled 
variable,  we  have  a  nonlinear  compensation  which  can  be  mechanized  in  a  relatively 
simple  way.  These  means  do  not  ordinarily  involve  severe  stability  problems.* 

In  the  weakly  coupled  case,  the  success  of  a  design  depends  heavily  upon  our 
knowledge  concerning  the  relationships  between  stability  derivatives  and  the 
parameters  which  change  them. 

After  sketching  these  two  different  approaches  (adaptive  control  and  non¬ 
linear  compensation)  let  us  return  to  our  concept  of  the  human  pilot  and  see  how 
the  principles  outlined  earlier  are  involved  in  the  choice  of  a  design  approach. 

The  most  significant  aspect  of  the  human  operator  is  his  adaptability. 

This  property  enables  him  to  control  a  system  with  widely  varying  parameters. 

Our  problem  with  the  peripheral  jet  ground  effect  machine  is  that  the  range 
of  stability  derivative  variations  exceeds  the  range  of  human  adaptability. 

If  we  knew  nothing  about  the  variations  of  the  stability  derivatives,  our 
only  choice  would  be  adaptive  control.  If  we  knew  everything  about  them,  non¬ 
linear  compensation  would  solve  the  problem  perfectly. .  We  are,  however,  in  a 
state  between  these  two  extremes.  The  actual  criterion  for  the  choice  can  be 
presented  in  the  following  question: 

Do  we  have  sufficient  information  about  the  variation  of  the  stability 
derivatives  so  that  the  use  of  nonlinear  compensation  and  feedback  results  in 
a  system  not  exceeding  the  range  of  adaptability  of  the  pilot  for  satisfactory 
control? 


*Stability  problems  may  arise  if  a  parameter  is  dynamically  coupled  to  the  con¬ 
trolled  variable  so  strongly  that  an  adjustment  of  the  feedback  loop  based  on 
this  parameter  must  be  considered  as  a  nonlinear  feedback  rather  than  as  a 
compensation. 
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The  most  important  point  in  this  question  is  that  we  wish  to  make  use  of 
the  adaptability  of  the  human  operator.  If  the  answer  to  the  above  question  is 
"yes",  then  there  is  no  need  to  go  into  the  design  of  a  more  complex  adaptive 
feedback  system.  The  answer  is  "yes."  Material  presented  in  this  report  pro¬ 
vides  sufficient  information  for  the  preliminary  design  of  a  nonlinear  compensa¬ 
tion  for  attitude  stability  augmentation  by  means  of  feedback  control. 

The  basic  design  principle  is  now  the  following:  We  endeavor  to  design 
a  feedback  loop  to  be  adjusted  at  any  height  in  such  a  way  that  the  pilot  is 
always  confronted  with  a  vehicle  which  is  easy  to  control.  We  cannot  achieve 
this  goal  perfectly  for  three  reasons: 

a.  We  cannot  predict  very  accurately  the  magnitudes  of  the  stability 
derivatives. 

b.  We  do  not  know  exactly  what  an  optimum  GEM  configuration  would  be  like 
from  the  point  of  view  of  handling  qualities . 

c .  We  have  made  numerous  approximations  in  order  to  simplify  the  system  to 
a  form  suitable  for  mathematical  synthesis. 

These  three  points  may  seem  discouraging  but  we  have  two  strong  points  to  make 
in  order  to  balance  all  the  uncertainties.  * 

a.  The  range  of  adaptability  of  the  pilot  can  absorb  the  uncertainties  of 
points  (a)  and  (c)  above. 

b.  Since  we  cannot  expect  to  hit  the  center  of  the  test  pilot's  range  of 
adaptability,  provisions  may  be  made  for  adjustments  of  the  feedback  gains  during 
test  flights.  Results  from  these  test  flights  can  then  easily  be  incorporated 
into  the  automatic  compensation  scheme. 

5.2.  Preliminary  design  of  a  feedback  system  for  attitude  stability  augmentation. 

Up  to  this  point,  our  effort  has  been  concentrated  on  simplifying  the  system 
while  keeping  the  most  significant  stability  problems  in  the  simplified  model, 


and  in  establishing  reasonable  and  practical  design  principles  for  the  feed¬ 
back  system. 

According  to  our  previous  discussion  we  are  not  going  to  consider  any 
lead  or  lag  contributed  by  the  pilot,  but  rather  leave  his  ability  to  provide 
such  contributions  as  a  safety  margin  which  can  be  expected  to  absorb  inaccu¬ 
racies  in  the  feedback  synthesis.  It  is  important  to  keep  in  mind  that  the 
primary  goal  of  the  pilot  is  to  steady  the  attitude.  The  use  of  attitude 
control  for  maneuvering  is  not  very  efficient.  Therefore,  the  task  of  fight¬ 
ing  disturbance  moments  is  strongly  emphasized.  Such  moments  are  caused  by  any 
change  of  the  center  of  pressure  under  the  base.  Experience  indicates  that 
small  changes  of  a  random  nature  always  occur.  At  heights  where  the  machine 
is  statically  unstable,  this  is  sufficient  to  initiate  an  exponential  divergence. 
Additional  disturbance  moments  are  caused  by  gusts . 

At  this  point  we  can  consider  three  basic  feedback  configurations.  The 
configuration  to  be  found  most  advantageous  will  be  investigated  in  further 
detail.  The  three  configurations  to  be  considered  are  the  following: 

a.  a  lead-lag  circuit  to  be  inserted  between  the  pilot  and  the  vehicle, 
ho  internal  control  loop  (see  Figure  26a); 

b.  internal  loop  with  rate  feedback  only  (see  Figure  26b)> 

c.  Internal  loop  with  combined  attitude  and  rate  feedback  (see  Figure  2 6C). 

Configuration  A  is  actually  an  open  loop  compensation,  and  the  response  to 

a  disturbance  moment  develops  uninhibited.  Among  the  three  configurations  this 
makes  the  task  of  the  pilot  the  most  difficult. 

In  configuration  B  the  effect  of  the  rate  feedback  is  to  slow  down  the  di¬ 
vergence  .  With  increased  damping  provided  by  the  rate  feedback,  the  unstable 

* 

root  moves  toward  the  origin.  However,  even  infinite  inner  loop  gain  cannot  move 
this  root  into  the  region  of  stability  for  unattended  operation. 
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It  should  be  mentioned  that,  if  there  were  no  additional  lag  in  the 
control,  we  could  count  on  the  pilot  to  provide  lead  sufficient  to  pull  the 
right  hand  side  root  over  into  the  region  of  stability.  However,  as  pointed 
out  earlier,  we  must  expect  a  lag  in  the  development  of  a  part  of  the  control 
moment.  Considering  this  additional  effect,  we  would  be  left  with  only  a 
narrow  range  of  pilot  gain  for  stability,  if  any .  This  approach  would  be 
likely  to  lead  to  a  highly  critical  design  and  would  leave  no  margin  for 
absorbing  the  inaccuracies  to  be  expected  in  the  design  information. 

Configuration  C  is  more  advantageous.  The  zero  of  the  transfer  function 
of  the  inner  loop  can  be  placed  on  the  left  hand  side  of  the  real  axis  with¬ 
out  any  contribution  from  the  pilot.  This  is  the  combined  effect  of  attitude 
and  rate  feedback.  Then  the  vehicle  which  confronts  the  pilot  can  be  made 
stable  in  attitude.  The  nature  of  the  control  which  he  exercises  changes 
from  moment  control  to  attitude  control.  Using  this  approach,  we  count  on 
the  pilot's  adaptability  only  in  order  to  counteract  the  effects  of  the 
presently  unknown  additional  lag  and  inaccuracies  in  the  design  information. 

A  single  vertical  gyroscope  can  provide  attitude  signals  for  Doth  the 
pitch  and  the  roll  control  feedback  loops .  This  is  a  very  reasonable  price 
for  the  advantage  of  Configuration  C.  Therefore  this  is  chosen  as  the  sug¬ 
gested  basic  feedback  configuration  for  the  attitude  stability  augmentation 
of  the  ground  effect  machine . 

The  transfer  function  of  a  feedback  path  (Figure  26c )  which  consists  of 
an  attitude  gyro  and  a  rate  gyro,  can  be  expressed  as: 

Gf  =  A  p  (  \  d  S  Tp  )  where  ~Tp  = 

This  indicates  that  the  location  of  the  zero  (  S  =  -  )  can  be  changed  by 

i  p 

adjusting  only  the  rate  feedback  gain,  /\  ^  .  For  an  adjustment  of  the  loop 


gain,  either  Ac  can  be  changed,  or  else  the  attitude  feedback  gain,  A  F  , 
and  the  rate  feedback  gain,  Ar  ,  must  be  changed  in  the  same  proportion. 
One  objective  of  the  design  is  that  the  product,  Ac  Act  be  as  large  as 
possible.  This  is  because  the  response  to  disturbance  moments  is  diminished 
almost  in  proportion,  without  any  contribution  from  the  pilot. 

It  has  been  shown  previously  that  the  left  hand  side  pole,  S|  ,  of  the 
attitude  transfer  function,  Ge  ,  (Equation  83)  does  not  change  much 
throughout  the  useful  height  range.  This  suggests  that  the  feedback  time 
constant,  —  Ar  /af  ,  can  be  kept  constant  so  as  to  approximately 
cancel  S|  at  all  heights:  *Tp  =  —  ^ 

The  inner  loop  gain  and  the  control  system  bandwidth  depends  to  a  large 
extent  on  the  additional  control  lag  time  constant,  Tc  .  We  can  arrive  at 
an  estimate  of  the  inner  loop  gain  from  the  following  considerations.  In 
order  to  obtain  a  simple  picture,  let  us  assume  that  we  have  managed  to  cancel 
the  root,  ,  and,  therefore,  we  are  left  with  two  transfer  function  poles: 

the  unstable  right  hand  side  GEM  pole,  S2  >  and-  the  control  lag  pole 

3c  =  ~  '/tc  •  As  a  criterion  for  the  gain  of  the  simplified  inner  loop 
we  suggest  a  damping  ratio,  =  ^=="  ,  as  a  convenient  and  reasonable 

compromise.  Much  higher  damping  is  not  easy  to  use  because  near  critical 

damping  )  the  roots  move  relative  rapidly  with  any  change  in  gain. 

We  are  not  able  to  be  exact  in  specifying  the  gain.  A  much  lower  damping 
would  not  provide  a  sufficient  margin  of  stability.  Although,  according  to 
experience,  it  might  be  expected  that  the  pilot  could  probably  handle  a  smaller 
margin  of  stability,  we  wish  to  stress  again  that  we  do  not  suggest  basing  the 
design  criteria  on  the  limits  of  the  pilot's  capability. 

Using  our  simplifying  assumptions,  and  writing  Ac/^\  +  s~Tc^  — 
we  arrive  at  the  following  block  diagram  for  the  inner  loop: 
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The  characteristic  equation  of  the  closed  loop  is 

s'2-  (S2t  sc)s  -t  Si  Sc  -V 


Aq  Af  Se  _ 

IS,0 


Since  the  roots  corresponding  to  a  damping  of  must  be  on  the 

lines  enclosing  an  angle  of  45°  with  the  negative  real  axis,  we  can  find  the 
inner  loop  gain  for  this  damping  by  substituting  S-  —  <3"-+  ^.00=.  —  CT-t-^CT 
in  the  characteristic  equation.  We  eliminate  C  from  the  two  equations  ob¬ 
tained  for  the  real  and  imaginary  parts,  both  of  which  must  become  zero  for 
the  same  (T  .  We  then  obtain: 


AjAf  =  S'SC  ['  +  (-fjO 


A  measure  of  the  bandwidth  of  the  inner  loop  is  \c5"\=  - 


(86) 


Sx  +■  Sc  1 


(rad/sec ) , 


These  results  indicate  how  strongly  the  control  lag  pole,  5c  >  influences  the 
adjustment  of  the  inner  loop  gain.  Examining  Equation  86,  we  find  that  if 

(Sc  1>  A  lS2|ma.*  }  the  influence  of  the  change  of  Sz  (i-e->  of  the  moment 
derivative)  on  the  desired  inner  loop  gain  is  fairly  small  (less  than  25$). 

This  means  that  if  the  ratio  of  the  control  gain  and  the  control  lag  root 
_  Ac/Scj=  Ac~Tc  ,  were  to  remain  constant  throughout  the  height  range,  our 


choice  of  feedback  configuration  and  dynamic  performance  criteria  would  probably 


satisfy  all  practical  purposes  with  a  fixed  feedback  gain,  A  F"  >  and  a  fixed 
rate  feedback,  Ar  .  However,  the  fact  that  the  desired  feedback  gain,  A  P" , 
is  inversely  proportional  to  Ac  ~Tc.  indicates  that  we  should  provide  a  variable 
gain. 

Since  the  variations  of  the  control  parameters  influence  the  controllability 
even  if  no  feedback  is  used,  it  would  seem  advantageous  to  provide  for  changes  so 
that  Ac  would  remain  nearly  constant.  That  is  to  say,  the  moment  versus  stick 
displacement  relationship  should  not  vary  appreciably  throughout  the  height  range. 
Considering  the  very  close  relationship  between  height  and  throttle  setting  for 
a  particular  configuration  it  seems  feasible  that  this  "automatic"  gain' adjust¬ 
ment  could  be  provided  purely  mechanically  by  changing  a  linkage-arm  ratio  by 
means  of  the  throttle.  It  is  then  conceivable  that  fixed  feedback  gains  could 
be  used  to  provide  satisfactory  dynamic  characteristics.  Another  advantage  of 
this  approach  is  that  the  gain  adjustment  would  be  an  open  loop  adjustment . 

This  eliminates  the  problems  of  measuring  height  and  obviates  any  problem  of 
dynamic  coupling  between  the  attitude  and  heave  modes  through  the  automatic 
control.  It  may,  however,  be  necessary  to  provide  for  some  adjustment  in  the 
gain  control  to  compensate  for  large  variations  of  the  load. 

It  may  also  be  pointed  out  that,  if  necessary,  an  additional  lead- lag  com¬ 
pensation  could  be  used  in  the  inner  loop  in  order  to  increase  its  bandwidth 
and  make  it  "tighter"  by  applying  higher  loop  gain.  Unfortunately  this  is  not 
immediately  possible  in  the  absence  of  information  available  about  the  vari¬ 
ation  of  the  control  system  parameters,  /Ac  and  ~Tc  .  In  the  present  state 
of  the  art,  therefore,  "tuning"  of  an  experimental  autopilot  throughout  the 
height  range  in  hovering  is  necessary  in  order  to  find  out  about  the  effect  of 
the  height  dependence  of  the  control  parameters.  Since  the  proposed  inner  loop 
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provides  stability  in  attitude  by  itself,  such  tests  can  easily  be  carried 
out  by  checking  moment-pulse  responses.  These  can  be  measured  by  recording 
the  attitude  gyro  output. 

The  previous  discussion  has  been  based  on  several  assumptions  and  simpli¬ 
fications.  It  is  appropriate  to  discuss  the  effects  of  inaccuracies  in  these 
assumptions  on  our  conclusions. 

In  our  simplified  approach  we  considered  three  real  roots:  the  two 
roots  S,  and  S2.  of  the  GEM  attitude  freedom  and  a  control  lag  root,  • 

We  have  arbitrarily  chosen  S|  to  be  approximately  cancelled  by  the  feedback 
zero  because  this  choice  offered  the  possibility  of  a  fixed  rate  gyro  gain. 

If  Sc.  should  turn  out  to  be  nearly  constant  throughout  the  height  range,  it 
could  be  cancelled  instead.  It  is  definitely  advantageous,  however,  to  cancel 
the  root  nearest  the  origin.  It  should  be  kept  in  mind  that  the  GEM  roots  Sf 
and  would  be  much  closer  to  the  origin  in  the  case  of  any  practical 

full  scale  design  than  they  are  for  the  eight  foot  model.  The  relative  radius 


of  gyratix  •>  of  a  practical  machine  would  be  larger  than  that  of  our  model.  In 
our  case  =  2.5  (  Pcj-  =  radius  of  gyration),  whereas  in  a  practical 

design  r/r>  <  2.0  seems  more  realistic.  With  the  value  =  2.0  and  a 

diameter  of  40  feet  instead  of  eight  feet,  the  GEM  attitude  roots  would  be  2.8 
times  closer  to  the  origin  than  the  roots  our  test  model  (if  the  relative 
nozzle  area  is  kept  the  same).  It  would  be  difficult  to  decide  whether  Si 
or  5c  would  be  closer  to  the  origin.  The  control  lag  can  be  expected  to  in¬ 
crease  with  size  but  to  decrease  i  1'  the  total  head  pressure  is  increased. 


The  following  practical  procedure  for  optimizing  the  location  of  the  feed¬ 
back  zone  can  be  suggested.  As  a  first  approximation  locate  the  zero,  (i.e., 
adjust  the  rate  gyro  gain)  with  the  assumption  that  the  relative  damping  ratio 
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does  not  change  with  size,  i.e.,  use  a  root  locus  plot  similar  to  Figure  23 
but  with  the  scale  changed  according  to  the  size  and  the  r<j.  ratio. 

Adjust  the  feedback  gain,  /Ap  ,  according  to  Equation  86  and  check  the  pulse 
response  at  several  heights.  Readjust,  Af  ,  for  approximately  critical 
damping.  Now  change  the  rate  feedback  (i.e.,  the  zero  location)  and  check 
for  uniformity  of  the  response  over  the  height  range .  Check  also  whether  a 
higher  feedback  gain  can  be  applied  for  the  same  damping .  Change  the  rate 
feedback  in  the  other  direction  and  repeat.  A  relatively  short  test  should 
lead  to  an  optimum  location  of  the  feedback  zero.  Now  keeping  this  zero  con¬ 
stant,  one  can  establish  experimentally  the  way  in  which  the  feedback  gain 
should  be  changed  in  order  to  compensate  for  the  change  in  the  control  pa¬ 
rameters  by  adjusting  for  approximately  the  same  damping  at  all  heights.  Since 
the  height  is  very  closely  related  to  the  throttle  setting  and  since  high  accu¬ 
racy  is  not  a  requirement,  the  gain  change  could  be  most  easily  instrumented  by 
gearing  appropriate  nonlinear  gain  potentiometers  for  the  pitch  and  roll  feed¬ 
back  to  the  throttle. 

We  should  like  to  emphasize  again  the  important  role  of  the  control  lag. 

We  can  conclude  that,  if  this  lag  can  be  avoided,  a  significantly  tighter  inner 
loop  with  a  considerably  larger  bandwidth  could  be  obtained  or,  if  such  improve¬ 
ment  is  not  needed,  the  elimination  of  the  attitude  feedback  and  the  application 
of  rate  feedback  alone  could  be  considered.  It  may  be  that  the  vehicle  designer 
can  exercise  some  control  over  this  feature  of  the  machine . 

We  have,  in  addition,  neglected  some  singulatities  of  the  inner  loop. 

We  assumed  that  the  servo  motor  transfer  function  poles  are  far  enough  removed 
so  that  their  effect  on  the  dynamics  of  the  controlled  vehicle  in  the  critical 
area  is  negligible.  We  have  shown  that  for  practical  machines  the  dominant  poles 
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are  much  closer  to  the  origin  than  those  of  our  eight  foot  test  model,  so  that 
this  assumption  is  certainly  justified.  We  have  further  neglected  any  effect 
of  dead  zone  or  backlash  in  the  controls.  This  would  be  very  difficult  to  con¬ 
sider,  unless  a  specific  design  was  analyzed.  We  have  also,  however,  neglected 
one  somewhat  favorable  singularity.  This  is  the  zero  which  appears  because, 
presumably,  only  part  of  the  control  moment  is  lagging .  If  we  write  the  con¬ 
trol  transfer  function  in  the  form 


-  A  O  t  1- 


Acz 

l  -V  S‘Tc 


Act -t  Ac2. t  sA  Cl  Tc  x 

\  +  S  Tc 


\  -V  S  (  Aci/Ac)To 

\  -v  s  Tc 


where  /Ac  ®AC!  /A  C.Z.  is  the  total  static  control  gain,  we  can  see  that 
the  control  zero  is  Ac/Ac.  times  as  far  from  the  origin  as  the  control 
pole.  It  can  be  expected  that  the  lagging  part  will  be  larger  if  the  largest 
part  of  the  weight  is  supported  by  the  base  pressure  and  pressure  differentials 
are  responsible  for  the  largest  part  of  the  moment.  The  control  zero  can, 
therefore,  be  expected  to  come  closer  to  the  control  pole  with  increases  in 
height,  but  it  could  reach  it  only  if  the  full  control  moment  was  to  develop 
without  delay.  Only  near  the  upper  limit  of  the  useful  height  range  could 
favorable  effect  of  this  zero  be  expected  to  be  noticeable.  This  should  actually 
have  no  influence  on  the  design. 

We  have  considered  the  control  of  the  GEM  in  the  vicinity  of  zero  attitude 
angle,  and  we  have  shown  that,  in  this  case,  the  coupling  between  the  attitude 
and  the  heave  degrees  of  freedom  can  be  neglected.  These  assumptions  are  quite 
valid  for  hovering.  Up  to  now,  we  have  discussed  the  synthesis  of  the  control 
system  for  hovering  only;  we  now  have  to  discuss  the  performance  we  can  expect 
of  our  system  in  forward  flight,  together  with  changes  or  additions,  which  seem 
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necessary.  Since  there  is  little  information  available,  at  the  present  time, 
concerning  the  stability  derivatives  in  forward  flight,  we  shall  not  be  able 
to  make  quantitative  predictions  as  in  hovering.  We  wish  to  emphasize  that 
we  are  concerned  here  only  with  the  control  of  a  GEM  at  low  speeds .  The  use 
of  more  or  less  conventional  aerodynamic  surfaces  for  control  is  not  con¬ 
sidered  because  of  the  wide  variety  of  possibilities  and  because  of  the  ex¬ 
istence  of  already  well  established  methods  for  the  analysis  of  such 
situations . 

Let  us  first  consider  forward  flight  with  the  machine  in  a  level  attitude. 
This  situation  has  been  discussed  in  the  section  concerning  forward  flight. 

From  the  point  of  view  of  our  control  system  we  have  to  consider  two  distinct 
problems . 

The  first  problem  arises  because,  as  the  velocity  is  increased,  an  in¬ 
creasing  moment  is  acting  on  the  GEM.  Rather  than  to  put  the  burden  of  counter¬ 
acting  this  moment  on  the  control  system  described  so  far,  it  is  advantageous  to 
counteract  at  least  a  large  part  of  this  moment  independently  of  the  automatic 
control  system.  This  could  be  achieved,  for  example,  by  allowing  for  a  moment 
arm  for  the  auxiliary  thrust  vector,  or,  alternatively,  by  gradually  shutting 
off  the  jet  at  the  front  edge  with  increasing  velocity.  The  latter  solution 
seems  more  desirable  because  it  tends  to  eliminate  the  moment  rather  than  fight¬ 
ing  it.  Also,  some  additional  lift  may  be  gained  by  adding  some  of  the  mass 
flow  saved  in  the  front  to  the  rear  edge .  The  disadvantage  of  this  approach 
is  the  more  complicated  control  system. 

After  the  problem  of  trimming  the  steady  state  moment  has  been  solved,  we 
still  face  a  second  problem:  the  changes  in  stability  and  control  derivatives 
with  velocity.  Because  of  the  scarcity  of  quantitative  information,  an  empirical 
approach  with  an  experimental  autopilot  is  suggested.  This  is  feasible  with  the 
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proposed  feedback  system  since  pulse  responses  can  be  used  for  the  tests. 

Once  the  desired  static  -control  gain  versus  velocity  characteristics  are  es¬ 
tablished,  we  can  use  an  approach  to  the  mechanization  similar  to  the  one 
already  employed  for  the  changes  in  height.  The  steady  state  velocity  is 
closely  related  to  the  throttle  setting  of  the  auxiliary  thrust  source,  there¬ 
fore  as  a  first  approach,  the  gain  adjusting  potentiometers  could  be  simply 
geared  to  this  throttle .  However,  it  can  be  expected  that  for  fast  accelera¬ 
tions  and  decelerations  the  throttle  setting  could  lead  the  velocity  too  much 
to  provide  smooth  control  during  the  acceleration.  If  this  is  the  case  a  lag 
should  be  inserted  between  the  throttle  and  the  control  gain  potentiometers. 

Concerning  the  effect  on  the  controlled  system  of  the  coupling  between  the 
attitude  and  the  translational  freedom,  the  following  argument  can  be  made. 
Since  we  have  a  rather  tight  attitude  control  only  small  deviations  from  the 
zero  attitude  angle  are  to  be  expected,  and  corrections  in  attitude  will  occur 
fast  enough  so  that  no  appreciable  change  in  velocity  can  occur.  In  other 
words,  the  forces  coupling  the  translational  freedom  to  the  attitude  mode  are 
kept  small  by  the  feedback  control.  We  can  say  the  same  thing,  however,  about 
the  moments  occurring  due  to  velocity  changes  and  gusts;  these  must  be  con¬ 
sidered  as  disturbance  inputs  to  the  attitude  control.  The  coupling,  however, 
can  be  considered  to  be  practically  unilateral,  and  no  additional  stability 
problems  are  to  be  expected. 

There  are  two  practical  situations  in  which  deviation  from  zero  attutide 
angle  must  be  desired..  One  is  banking  in  a  turn,  the  other  is  enhancing  de¬ 
celeration  and  acceleration  by  tilting  the  vehicle  in  pitch.  Since  we  assume 
that  the  height  at  which  the  GEM  flies  is  determined  by  the  obstacle  to  be 
cleared,  tilting  can  be  used  only  if  reserve  power  is  available  so  that  the 
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lowest  edge  remains  at  the  reference  height.  It  is  obvious  that  only  small 
angles  can  be  achieved  in  practice.  Because  of  the  attitude  feedback  the  equi¬ 
librium  angle  can  be  commanded  by  the  control  stick  input .  It  may  be  expected 
that,  for  the  practical  small  angles,  the  stability  derivatives  do  not  change 
appreciably.  The  coupling  between  the  heave  and  the  attitude  mode  is  small 
with  small  angles,  and,  since  the  heave  mode  is  inherently  stable  and  damped, 
no  severe  stability  problems  are  expected  in  connection  with  small  tilt  angles. 

In  summary,  the  relatively  simple  autopilot  devised  on  the  basis  of  a 
fairly  detailed  analysis  of  the  hovering  case  can  be  expected  to  also  work 
satisfactorily  in  forward  flight,  after  it  has  been  provided  with  suitable 
automatic  gain  adjustments.  In  the  present  state  of  the  art,  the  gain  adjust¬ 
ment  functions  have  to  be  established  experimentally.  Once  these  functions 
have  been  established  they  can  be  permanently  set,  and  there  is  apparently  no 
need  for  adaptive  control.  Great  precision  in  the  synthesis  of  closed  loop 
system  dynamics  is  not  required  since  we  can  count  on  the  pilot's  ability  to 
handle  a  considerable  range  of  vehicle  characteristics. 


6.  CONCLUSIONS 


As  is  usually  the  case  when  one  attempts  to  synthesize  a  feedback 
control  system,  the  large  part  of  the  effort  is  this  investigation  went 
into  the  analysis  of  the  system  to  be  controlled.  In  the  course  of  an 
effort  to  isolate  the  most  significant  problems,  the  stabilization  of 
the  ground  effect  machine  in  attitude  was  singled  out.  In  the  theoreti¬ 
cal  approach  the  principle  of  the  balance  of  momentum  changes  and  forces 
was  used.  Derivations  of  expressions  for  the  undamped  natural  frequency 
and  the  damping  in  heave  as  well  as  for  the  static  attitude  moment  yielded 
results  which  reveal  the  basic  mechanisms  influencing  these  parameters  and 
which  are  in  good  agreement  with  experimental  results.  This  agreement  es¬ 
tablished  a  certain  amount  of  confidence  in  the  theoretical  results  so 
that  they  could  be  used  in  the  preliminary  design  of  a  feedback  control 
system.  A  discussion  indicated  the  influence  of  different  physical  pa¬ 
rameters  as  well  as  the  influence  of  scaling  on  the  theoretical  results . 

Especially  interesting  is  the  form  of  the  result  for  the  attitude 
moment  derivative  which  indicates  the  sources  of  both  the  stabilizing 
and  the  destabilizing  static  moments.  Both  of  these  moments  can  be 
attributed  to  losses  in  the  base  area:  the  stabilizing  moment  to  a 
loss  in  base  pressure  (due  to  vortices  and  other  possible  causes),  the 
destabilizing  moment  to  the  loss  of  momentum  of  the  crossflow  in  the  base 
area  (primarily  due  to  dissipation).  The  change  of  the  balance  of  these 
moments  against  each  other  causes  the  crossover  from  static  stability  to 
instability.  This  result  has  been  obtained  by  way  of  some  rather  crude 
approximations;  nevertheless,  it  seems  to  provide  a  frame  into  which 
more  refined  analysis  of  the  future  can  be  expected  to  fit. 
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Experiments  were  carried  out  with  an  eight  foot  diameter  model  to 
yield  information  about  the  static  and  the  dynamic  stability  derivatives 
in  a  hover  condition.  In  general,  satisfactory  agreement  with  the  theo¬ 
retical  results  was  obtained.  However,  from  the  measurement  of  attitude 
damping,  unexpected  results  were  obtained  for  which  no  theoretical  expla¬ 
nation  could  be  given.  These  tests  indicated  the  presence  of  a  degree  of 
damping  sufficiently  large  to  influence  considerably  the  preliminary  de¬ 
sign  procedure  of  the  feedback  system.  This  influence  seems  significant 
enough  to  suggest  exploration  of  the  effect  of  scaling  on  attitude  damping. 
It  was  not  explored  in  the  present  study. 

In  the  course  of  determining  design  principles  for  the  automatic  sta¬ 
bility  augmentation,  it  was  decided  that  the  adaptability  of  the  human 
pilot  could  be  effectively  employed.  This,  and  the  consideration  of  dis¬ 
turbance  moments,  leads  to  the  selection  of  a  combined  attitude  and  rate 
feedback  control  system  with  open  loop  gain  adjustments  to  take  care  of 
the  significant  effects  of  flight  condition  variations. 

The  preliminary  synthesis  brings  out  the  significance  of  the  moment 
control  lag  and  the  variation  of  the  moment  control  effectiveness  with 
height.  Because  of  the  significant  role  of  these  parameters  the  study  of 
the  attitude  control  dynamics  of  practical  GEM  control  systems  can  be  sug¬ 
gested  as  a  worthwhile  subject  for  future  investigations.  A  series  of 
arguments  leads  to  the  conclusion  that  the  suggested  feedback  control  con¬ 
figuration  could  be  expected  to  perform  satisfactorily  in  hovering  and  in 
forward  flight,  up  to  speeds  where  aerodynamic  control  surfaces  become 
effective.  This  is  the  case  even  though  very  little  information  is  availa¬ 
ble,  at  this  time,  concerning  the  variation  of  stability  parameters  in 
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forward  flight.  Therefore,  the  installation  of  an  experimental  auto¬ 
pilot  built  according  to  the  preliminary  design  described  in  this  re¬ 
port,  and  installed  in  a  piloted  flying  model,  is  suggested  as  the 
next  important  step  in  the  study  of  automatic  stability  augmentation. 

It  should  be  kept  in  mind  that  this  report  is  concerned  only  with 
the  overland  characteristics  of  the  ground  effect  machine.  Experimental 
over-water  studies  have  indicated  a  different  behavior  of  the  stability 
derivatives.  The  suggested  feedback  control  configuration,  however,  is 
able  to  cover  such  a  wide  range  of  characteristics  that  a  different  gain 
adjustment  program  can  be  expected  to  provide  satisfactory  stability 
augmentation  over-water  also.  Switching  from  one  set  to  another  set  of 
gain  adjusting  potentiometers  on  the  part  of  the  pilot  would  probably  be 
preferred  to  an  automatic  adaptive  control  system  because  of  the  faster 
adaptation.  The  application  of  adaptive  control  principles  might  be 
needed  to  replace  the  human  adaptivity  if  remote  control  of  ground  effect 
machines  were  contemplated. 
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APPENDIX  A 


EFFECT  OF  GYROSCOPIC  COUPLING 

In  order  to  be  able  to  estimate  the  effect  of  gyroscopic  coupling, 
let  us  investigate  the  following  case.  Let  us  assume  that  a  sinusoidal 
oscillation  takes  place  around  the  roll  axis,  with  some  frequency  OOxo  , 
which  cannot  be  easily  controlled  by  the  pilot.  Let  the  amplitude  of 
this  oscillation  be  /\x  •  A  sinusoidal  gyroscopic  moment  around  the 
pitch  axis  induces  an  oscillation  of  the  same  frequency  with  an  ampli¬ 
tude  V  The  amplitude  ratio  A*/a,  can  be  considered  as  a  practi¬ 
cal  measure  of  the  significance  of  the  gyroscopic  coupling.  It  can  be 
easily  shown  that 

A if  1*  -ax  (A-i) 

A*  Tif- 


where  is  the  moment  of  inertia  of  the  propeller  and  shaft  aligned 

with  the  Z»axis,  ru  is  the  angular  velocity  of  the  propeller  and  Ti^. 
is  the  moment  of  inertia  of  the  machine  around  the  axis .  If  we 
choose  C =  IfT/sec  as  the  "input"  frequency.  Equation  A-l  can  be 
written  as 


where 


is  the  propeller  speed  in  revolutions  per  second. 


(A-l-a) 
If  we  sub¬ 


stitute  in  this  expression  the  measured  values  of  the  tested  eight  foot  model, 
'3x  =  .07  (ft.  lb.  sec.2);  ~Jxj-  =  13*4  (ft.  lb.  sec.2)  we  find  with  the  max¬ 
imum  propeller  speed(4-350  RFM)  A*f  /  Ap  =  O.38.  This  amplitude  ratio  is  not 
negligible . 
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The  following  discussion  indicates  the  effect  of  scaling  on  the  gyroscopic 
coupling.  For  two  machines  of  similar  shape  and  weight  distribution,  but  differ¬ 
ent  in  size,  the  moments  of  inertia  differ  by  an  approximate  ratio 


3m-z  VA_d/  Pl-x  D*  (a-2) 

3^i  v/|  D,1  Pt-i 

where  R_i  and  are  "the  specific  loads  in  lb/ft  .^.  If  the  index  2 

marks  the  larger  machine,  the  moment  of  inertia  ratio  increases  by  more  than 
the  fourth  power  of  the  diameter  ratio  if  the  larger  machine  also  has  a  larger 
specific  loac  . 

In  order  to  make  an  estimate  of  the  increase  in  the  moment  of  inertia  of 
the  propeller,  the  horsepower  requirement  must  be  considered.  The  horsepower 
required  is  proportional  (  indicated  by  )  to  the  product  of  the  mass  flow 
and  the  square  of  the  jet  velocity. 


3/-. 

UP  rv  ~  m-Tor  VTy  ^  SV»  P T  ^A"3^ 

where  is  the  nozzle  area,  P-p  the  total  pressure,  M" 3-  is  the  jet 

velocity  and  is  the  total  mass  flow  of  the  jet.  By  definition  the 

lift  is  L=V7-  ATtot  (where /A  is  the  augmentation  factor  and  3tcu  is 
the  total  jet  momentum).  From  this  expression  we  obtain  Pu  S  'w  A  5m  PV 

by  using  Pu  =  &nd  fac-t  that  T  ^  SA  Pt  •  We  can  use  this 

o 

relationship  to  eliminate  P-p  from  the  horsepower  relationship  giving, 


(A-M 
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The  propeller  horsepower  can  be  scaled  approximately  as  HP*->  Rprtrp 
or  as 

HPpv,  ~  Rptp  U-5) 

if  the  tip  speed  is  unchanged.  Comparing  the  last  two  relationships  we 
find  that  the  propeller  diameter  should  be  proportional  to  the  GEM  di¬ 
ameter  in  the  case  of  a  circular  machine,  if  the  scaling  is  made  for  the 
same  augmentation  and  specific  load.  We  assume  that  the  mass  of  the  pro¬ 
peller  increases  only  approximately  with  the  square  of  the  diameter  be¬ 
cause  relatively  less  structural  strength  is  needed  because  of  smaller 
RPM  of  the  larger  propeller.  Then 


Sf  ~  (4;)2  (£)*(%/ 


Applying  the  formulas  A-6  and  A-2  to  A-l-a 


(A-6) 


Using  this  formula  with  N**  -  I  the  gyroscopic  coupling  of  the 

Princeton  twenty  foot  model  could  be  characterized  from  the  data  of  the 
eight  foot  model  with 

0.5 

Actually,  however,  the  propeller  diameter  is  larger  by  a  factor  of  only  two 
and  the  necessary  horsepower  is  obtained  by  means  of  a  larger  tipspeed.  This 
leads  to  a  corrected  estimate  of  ="  0.18.  Considering  that  our 

fixed  frequency  CO*0  =  <XfT/sec.  is  rather  fast  for  an  attitude  oscillation. 


9k 


this 


,^2/^*2  ratio  indicates  a  comparatively  small  effect  of  the  gyro¬ 
scopic  coupling.  This  is  confirmed  by  the  pilot's  opinion.  According  to 
Equation  A-7  this  situation  can  be  expected  to  remain  unchanged  as  long  as 
the  increase  of  the  specific  load  is  not  larger  than  the  square  root  of  the 
increase  in  diameter. 

This  discussion  leads  to  the  following  general  conclusions: 

a.  The  significance  of  gyroscopic  coupling  in  practical  machines  can 
be  considered  practically  negligible  except  for  single  rotor  ground  effect 
machines  with  very  large  specific  loads. 

b.  The  scaling  of  the  gyroscopic  coupling  must  be  considered  if  the 
dynamics  of  a  single  propeller  ground  effect  machine  are  to  be  investigated 
on  a  scaled  model. 

c.  In  the  case  of  large  specific  loads,  gyroscopic  coupling  can  and 
should  be  practically  eliminated  by  using  a  number  of  contrarotating 
propellers . 
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APPENDIX  B 


EXPERIMENTAL  INVESTIGATION 

The  discussion  which  follows  describes  the  model  used  for  the  tests, 
the  testing  procedure  and  the  results  of  the  tests.  Because  of  limitations 
in  the  length  of  time  available,  no  tests  were  made  in  forward  flight. 

Model  and  Test  Setup:  The  model  on  which  the  experiments  were  carried 
out  is  a  circular,  annular  jet  ground  effect  machine  with  a  diameter  of  8  feet 
and  a  weight  of  170  pounds .  The  jet  is  one  inch  thick  and  is  inclined  inward 
at  an  angle  of  4  5  degrees.  Power  is  supplied  by  a  four-blade  propeller  driven 
by  a  three-phase,  115  volt  electric  motor.  The  propeller  is  2k  inches  in 
diameter.  Power  to  the  motor  is  controlled  by  means  of  a  Variac  located  at 
the  control  station. 

The  model  is  attached  to  a  rig,  as  shown  in  Figure  28,  so  that  it  has 
two  degrees  of  freedom.  The  bearings  supporting  the  model  are  located  on  an 
axis  passing  through  the  center  of  gravity.  This  allows  freedom  in  pitch 
about  this  axis  only.  The  bearings,  in  turn,  are  supported  by  the  yoke  shown 
in  Figure  27.  A  vertical  shaft  from  the  yoke  extends  through  a  sleeve  in  the 
main  rig.  This  allows  freedom  of  motion  in  the  heave  mode.  The  model  can  be 
locked  in  either  heave  or  pitch  by  means  of  notched  shafts  extending  from  the 
yoke  in  the  case  of  the  pitch  lock,  and  from  the  main  rig  in  the  case  of  the 
heave  lock  (see  Figure  27).  The  notches  are  made  at  intervals  of  one  degree 
in  pitch  and  one  inch  in  heave.  A  pneumatic  release  system  effects  the  re¬ 
lease  of  the  model  in  either  pitch  or  heave  or  both  with  power  on.  (Springs 
which  can  be  seen  in  the  Figure  pull  the  notched  shafts  free  when  released.) 

Heave  and  pitch  motion  are  read  by  means  of  potentiometers .  Referring 
to  Figure  27,  the  pitch  pot  may  be  seen,  with  its  linkage,  at  the  base  of 
the  yoke.  The  pot  is  rigidly  attached  to  the  yoke  while  its  wiper  arm  is  free 
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to  move  with  the  model.  The  heave  pot  may  be  seen  in  Figure  27  just  above 
the  yoke  at  the  apex  of  the  linkage  connecting  the  main  rig  and  the  yoke . 

This  pot  is  rigidly  attached  to  the  main  rig  while  its  wiper  arm  is  free  to 
move  with  the  yoke.  The  wiper  voltages  are  fed  into  a  Sanborn  recorder  at 
the  control  station  and  also  into  meters  for  direct  reading. 

The  weight  of  the  machine  can  be  increased  by  attaching  copper  weights 
to  the  structure  and  can  be  decreased  by  means  of  a  spring  attached  to  a 
hoist  above  the  model.  The  spring  is  sufficiently  soft  so  that  it  has  only 
a  small  effect  on  the  period  of  the  transient  oscillations  in  heave.  Cor¬ 
rections  are  made  to  account  for  the  fact  that  the  mass  does  not  decrease 
when  the  spring  is  used  to  reduce  the  weight. 

Static  moments  are  measured  by  means  of  a  strain  guage  which  is  bonded 
to  a  beam  attached  to  the  model  at  one  end  and  to  a  cable  which  is  suspended 
fromabove  at  the  other  end.  Weights  are  added  to  this  end  of  the  model  so 
that  all  moment  readings  are  positive.  The  signal  is  read  on  the  output  of 
a  Sanborn  amplifier  at  the  control  station. 

Testing  Procedure;  Measurements  were  made  to  show  the  variation  of  the 
lift  force  with  height  for  various  constant  values  of  the  motor  KPM.  Weight 
was  both  added  to  and  subtracted  from  the  basic  weigh  in  20  pound  increments. 
For  each  weight,  the  resulting  hovering  height  was  recorded  at  three  specific 
values  of  the  motor  RPM.  The  results  of  these  measurements  are  shown  in 
Figure  2.  The  curves  of  Figure  3  show  the  effect  of  small  attitude  changes 
on  the  lift . 

A  series  of  experiments  was  conducted  to  determine  the  dynamic  behavior 
of  the  model  in  heave.  To  accomplish  this,  the  model  was  locked  in  pitch  so 
that  it  had  only  vertical  freedom.  Because  of  the  presence  of  some  flexibility 
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in  the  supporting  structure  the  model  tended  to  exhibit  some  attitude  change 
with  variations  in  the  power  supplied  to  the  propeller;  however,  these  changes 
were  usually  small  and  were  accounted  for  whenever  possible. 

The  tests  were  conducted  as  follows:  The  model  was  raised  to  a  certain 
height,  the  attitude  locked  at  a  certain  angle  and  the  model  weight  set  by 
either  adding  weights  or  increasing  the  spring  force.  The  propeller  was  then 
brought  up  to  an  RPM  which  corresponded  to  an  equilibrium  height  either  some¬ 
what  greater  or  somewhat  less  than  the  one  at  which  the  model  was  set.  The 
model  was  then  released  and  the  resulting  transient  was  recorded.  The  fre¬ 
quency  and  damping  ratio  were  then  determined  from  the  recorded  trace.  System¬ 
atic  variations  were  made  in  the  model  weight,  pitch  angle,  equilibrium  height 
and  also  in  the  size  of  the  input.  It  was  found  that  inputs  of  about  plus  or 
minus  1.5  inches  were  a  good  compromise.  (Smaller  inputs  did  not  give  a  suf¬ 
ficient  number  of  readable  overshoots  for  determining  the  frequency  and  damping, 
and  larger  inputs  would  run  the  risk  of  causing  nonlinear  behavior).  The  re¬ 
sults  of  the  tests  are  shown  in  Figure  6  and  in  Figures  9  to  12. 

Pitch  Freedom  Tests;  The  desired  information  in  the  pitch  degree  of  free¬ 
dom  was  somewhat  more  difficult  to  obtain  because  of  the  fact  that  the  model 
was  unstable  in  pitch  at  all  but  the  very  lowest  heights  .  Essentially  two 
measurements  were  made  at  each  height.  The  variation  of  pitching  moment  with 
pitch  angle  was  recorded,  and  the  variation  of  pitch  damping  was  determined 
separately.  In  each  case  the  model  was  rigidly  fastened  at  a  given  height 
from  the  ground.  To  determine  the  effect  of  variation  of  the  machine  weight, 
each  test  was  made  at  three  different  values  of  the  motor  RPM.  This  was  equi¬ 
valent  to  varying  the  weight  of  the  machine  since  it  had  no  vertical  freedom. 


As  described,  previously,  the  moment  measurements  were  made  using  a  strain 
gauge  which  was  bonded  to  a  beam  at  the  edge  of  the  machine .  The  machine  was 
set  at  the  desired  height,  the  cable  adjusted  to  give  the  desired  pitch  angle 
and  then  the  propeller  was  brought  up  to  speed.  Moment  readings  were  taken 
at  each  of  three  previously  determined  RPM  readings.  The  results  of  these 
tests  are  shown  in  Figure  21.  In  all  of  the  moment  measurements  there  was 
a  considerable  amount  of  fluctuation  present  and  each  of  the  readings  repre¬ 
sents  an  average  taken  over  several  seconds  of  the  motion  of  the  needle  in¬ 
dicating  the  strain  gauge  output.  The  fluctuations  appeared  to  be  due  to  the 
unsteadiness  of  the  flow  beneath  the  model. 

Pitch  damping  was  determined  by  a  test  in  which  a  fixed  moment  was  applied 
to  the  model  at  different  height  and  RPM  settings.  The  time  history  of  the 
pitch  angle  was  recorded  until  the  model  struck  the  ground.  The  moment  was 
applied  in  such  a  manner  as  to  insure  that  the  pitch  angle  varied  through  the 
linear  portion  of  the  moment  curve.  Assuming  that  this  motion  could  be  des¬ 
cribed  by  a  second-order  differential  equation  in  which  the  inertia  and  spring 
force  (negative)  were  known,  curves  showing  the  pitch  angle  response  to  a  step 
input  as  a  function  of  damping  ratio  were  plotted.  From  these  curves  the  damp¬ 
ing  ratio  corresponding  to  the  measured  pitch  angle  could  be  found.  The  varia¬ 
tion  of  the  damping  ratio  with  height  is  shown  if  Figure  22. 
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FIGURE  II  EXPERIMENTAL  HEAVE  DAMPING 
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FIGURE  22  EXPERIMENTAL  PITCH  DAMPING 
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FIGURE  24  EFFECT  OF  RPM  ON  PITCHING  MODE  (EXPERIMENTAL) 


FIGURE  25  EFFECT  OF  DAMPING  ON  STATICALLY 
UNSTABLE  SECOND  ORDER  SYSTEM 


FIGURE  26  FEEDBACK  CONFIGURATIONS 


the  Model 


I 


AIART  PROGRAM 
Technical  Report 
Distribution  List 


ADDRESS 

1.  Chief  of  Transportation 
Department  of  the  Army 
Washington  25,  D.C. 

ATTN:  TCACR 

2 .  Commander 

Wright  Air  Development  Division 
Wright -Patters on  Air  Force  Base,  Ohio 
ATTN:  WCLJA 

* 

3 .  Commanding  Officer 

U.S.  Army  Transportation  Research  Command 
Fort  Eustis,  Virginia 
ATTN:  Research  Reference  Center 
ATTN:  Aviation  Directorate 

4.  U.S.  Army  Representative 
HQ  AFSC  (SCS- 3  ) 

Andrews  Air  Force  Base 
Washington  25,  D,C. 

5 .  Director 

Air  University  Library 
ATTN:  AUL-8680 

Maxwell  Air  Force  Base,  Alabama 

6.  Commanding  Officer 
David  Taylor  Model  Basin 
Aerodynamics  Laboratory 
Washington  7,  D.C. 

7-  Chief 

Bureau  of  Naval  Weapons 
Department  of  the  Navy 
Washington  25,  D.C. 

ATTN:  Airframe  Design  Division 
ATTN:  Aircraft  Division 
ATTN:  Research  Division 

8.  Chief  of  Naval  Research 
Code  46l 

Washington  25,  D.C. 

ATTN:  ALO 


NO.  OF  COPIES 


2 


4 

3 


1 


l 


l 


1 

l 

1 


1 


ADDRESS 


NO  =  OF  COPIES 


9-  Director  of  Defense  Research  and  Development 
Room  3E  -  IO65,  The  Pentagon 
Washington  25,  D.C. 

ATTN :  Technical  Library  1 

10.  U.S.  Army  Standardization  Group,  U.K. 

Box  65,  U.S.  Navy  100 

FPO  New  York,  New  York  1 

11.  National  Aeronautics  and  Space  Administration 
1520  H  Street,  N.  W. 

Washington  25,  D.C. 

ATTN :  Bertram  A .  Mulcahy 

Director  of  Technical  Information  5 

12 .  Librarian 

Langley  Research  Center 

National  Aeronautics  &  Space  Administration 

Langley  Field,  Virginia  1 

13 .  Ames  Research  Center 

National  Aeronautics  and  Space  Agency 
Moffett  Field,  California 

ATTN :  Library  1 

14.  Armed  Services  Technical  Information  Agency 
Arlington  Hall  Station 

Arlington  12 ,  Virginia  10 

15 •  Office  of  Chief  of  Research  and  Development 
Department  of  the  Army 
Washington  25,  DC. 

ATTN:  Mobility  Division  1 

16.  Senior  Standardization  P-^Dret  entative 
U.S.  Army  Standardization  Group,  Canada 
c/o  Director  of  Weapons  and  Development 
Army  Headquarters 

Ottawa,  Canada  1 

17.  Canadian  Liaison  Officer 
U.S.  Army  Transportation  School 
Fort  Eustis,  Virginia 


3 


ADDRESS 


NO.  OF  COPIES 


18.  British  Joint  Services  Mission 
(Amy  Staff) 

DA QMS  (Mov  &  Tn) 

1800  "K"  Street,  N.  W. 

Washington  6,  D.C. 

ATTN:  Lt  Col.  R.J.  Wade,  R.E.  3 

19.  Office  of  Technical  Services 
Acquisition  Section 
Department  of  Commerce 

Washington  25,  D.C.  2 

20.  Librarian 

Institute  of  the  Aeronautical  Sciences 
2  East  64-th  Street 

New  York  21,  New  York  2 

21.  Chief 

U.S.  Army  Research  and  Development  Group  (Europe) 

ATTN:  USATRECOM  Liaison  Officer 
AFO  757 

New  York,  New  York 


1 


_Ac  cess  ion  No.  _  UNCLASSIFIED  AD _ Accession  No. 


p 

o 

<D 

P  I 

P  >> 
W  co  P 
CD  P 
d  C  H 
2  P  P 

g  a  ,g 
O  o  as 
P  eg  P 
C2M 


P 

cj 

£ 

P 

CJ  I  P 


d  O  Id  C 

6  p  g  p 
o  -p  3  ^ 
p  a  o  o 
3  O  f(  3 
<<  o  o  s 


$  £ 
as  iH  -3- 

P  I  CM 
P  -4"  LCN 

a  -d  i 
o  <;  o 
O  Q  Eh 


fe  H  O 
O  l 


-p 

-P 

o 

O 

< V 

<u 

Cm 

1 

<H 

<M 

>> 

O 

1 

CM 

W 

w 

-p 

•H 

W 

(D 

•rH 

■P 

rH 

q 

rH 

<3 

O 

nd 

81 
ON 
CM  Sc 
c/N  Eh 

i 

o  *> 

EH  O 

I  O 

t—  o 

t-  I 

rl  rl  p 

i  o  T, 

-d  l  O 

4(0  ft 
•3;  co  cu 
oik 

ON 

•  d 

O  •  4) 

K  o  P 

a  p 

+3  P 

O  P  10 
as  o  co 
P  <D  aS 
■P  H-J  H 
d  O  a 


O  P  3  d 

P  C!  O  Cl 
3  O  P  3 
<002 


P  t- 
ci  C— 
aJ  H  -4 
P  I  CM 
P  -4-  u~\ 
d  -d  l 
O  <  U 
OPE 


P 

o 

<D 

Pi  I 

P  >, 
W  co  P 
4)  P 
3  rl 


P 
Cl 
Cl 
Pi 
I  Pi 
W  CO 

h  ai 

D  d  d 


P  rH 

I  P  C 
O  P  P 

p  a  o 
d  o  p 
<  o  o 


K  EH  O 

O  I 


<  o  o 


p  c— 

cl  t— 
aJ  rH  . 

P  i 

S3 


fe  Eh  O 

O  i 


fe  Eh  O 

O  I 


I  O  P 
-d  l  O 
-4-  CO  Pi 
cn  a> 
Q  I  K 
ON 

•  d 

o*ai 
a  o  p 
a  Pi 

p  P 
O  P  CO 
a!  ci  co 
P  a>  a) 
p  "~5  P 
POo 
O  P  d 
U  P  D 


Project  No.  9-38-01-000,  TK902  TC-524  Project  No.  9-38-01-000,TK902  TC-$2k 

Unclassified  Report  Unclassified  Report 


a  g 

O  ft 

•H  O  ft 
ft  ft 
aS  Ft 

N  •>  O 
■h  a  <m 
ft  O 
•H  ft  W 
ft  ft  G 
n)  cd  O  1 
-P  G  ft 
to  g  to  * 
Pt  to 
GOG 

ft  -ft  Pt 

G  §  [3  •! 


o  >  <L) 
G  O  ft 
P  P 

G  O  ft 
O  ft  O 

o  o 

*■£3  ft 
G  G 
CJ  P  > 


•0  t* 
>?  P  o 

t)  p  H 

0  to 
p  0 
M  G  G 

tdg  * 


•S  §  % 
S  o)  5 


0  - 
t  G  g  co 

)  o 3  G 

aJ  G 
i  to  cd 
G  P 

*  g  ft  rt 
P  60  o 
G  ft 


G  ft  G 
P  ft  P 
ft  cd  to 

>>  G 
G  a)  to  aS 


1  cd  P  <u 
o  o  ft 

Dft  ft  X 
CO  G 
>5  ft 

£  °  ft 
ft  O 
to 

.  ft  G  to 
O  O  0) 
•pH  irH 

tu  p  p 

O  <0  G 
G  0,  to 
01  ft 
G  4)  cd 


C_l  t_, 

P  O  P  0 

c  a  x: 

0  bp  O  P 
SCO 
0  *pH  ^ 
D  ft  ^  O 
P  £?  CJ 

bp  cd  cd  >3 


co  co 
"d  0 
i  co  0  bp 
•H  -p  G 
a  cd 


0  P 
•  0 

CO  0  -p 
CO  ,Q  0 


'  <HH  0  P 

CO 

0  P 

£0  p 
0  (0 
CO  0 
•  CQ  P 
0  P 
O  O  0 

rH  Tj 

cd 

ia  co 
•h  ^ 

^  M  Pc 

$53 

£ 

O  CJ  tu 
8  to  tl 


t  G  8  o 
MS  G  >j 
id  /i  P 
0  ft 
>5  ft  4)  H 
p  rH  4)  *H 
O  ft  ft  ft 
P  G  G 
O  P  4)  p 

G  4>  p  ft 
ft  d  i  G 
to  ft  G  0 
•H  CO  G 
ft  G  0 
G  O  G  4) 
co  O  G  ft 
p 

G  >  4) 

H  O  d  to 
G  3  G 

4)  tO  ft  ft  • 

P  ft  P  , 

G  O  p  4) 

G  P  0 

rH  ft  G  ft  - 

4)  G  G 
■ft  G  G  • 

O  O  r 

B  G  o 


O  ft  G  ft 


ft  ft  0 
ft  ft  P 


i — f  0  ft 
O  ‘ft  p 
P  ft  G 
p  G  G 


G  a 
p  too 
G  o 
<D  ft 

a  o  t>, 

OOP 

o  a  ft  o 

p 

4>  G  O 
G  4)  G 
ft  Pt  ft 

O  G 
ft  ft 

°  s3  Id 

G  G  G 
O  P 
ft  G 
P  ft  ft 
)  G  tU  ft 
ft  P  > 
P  G  O 
0  G  G 
>  bO  ft 


l 

ft  ft 
<U  ft  O 

G  ft 
ft  Ct3  G 


ft 

C  G  CO 
G  P 
G  G 


i  ft  a 

G  P  G 

P  ft  P 

ft  G  G 


N 

o 

A 

•H 

P 

bC 

l  o 

P  X 

3  B 

rC 

3 

•H 

P 

S' 

•H 

a 

cp 

-p 

O 

0 

•H 

•H  -f. 

5  H  p 

P 

H 

o 

3  ! 

r— 1 

o 

G 

§ 

rH 

CO 

> 

O  X5 

G 

P 

0 

G 

•H 

•H 

CO 

p 

•H 

0 

<P 

p 

4. 

i  P 

rH 

M 

pi 

o 

-P 

G 

o 

u 

P 

G 

P  C 

>  P  0 

CO 

o 

K 

cd 

cd 

O 

<HH 

Ph 

cd 

0 

G 

G  rG 

p 

0 

p 

-P 

u 

•H 

Ph 

P 

e 

0  b 

a  c 

OOP 

S 

p 

• 

0 

CO 

K 

CO 

0 

cd 

•H 

r  o 

P 

G 

CO 

0 

P 

CO 

CJ 

rG 

rH 

5 

P 

0  -H 

l  Cm 

cd 

O 

CO 

P 

0 

O 

0 

cd 

P 

cd 

p 

0 

G  O 

G 

o 

0 

B 

rC 

rH 

p 

O 

o 

Ph 

p  a 

O 

•H 

G 

G 

3  1 

pH 

G 

•H 

cp 

X 

sp! 

cd  >3 

6 

p 

0 

3 

p 

g 

X 

o 

CO 

0 

p  p 

•H 

G 

> 

o 

cd 

u 

cd 

W 

•H 

•H 

ft 

rH  *H 

H 

•H 

Ph 

o 

rH 

P 

CO 

ft 

O 

Cm 

:>>  0 

0  r—l 

0 

a 

P 

CO 

4h 

O 

G 

Pt 

O 

p  i — i 

0  *H 

P 

o 

CJ 

p 

P 

u 

• 

S 

ft 

CO 

O  rQ 

ft  ft  Oh 

a 

0 

G 

O 

rH 

0 

CO 

G 

co 

p  cd 

cd 

<M 

0 

cm 

O 

k 

0 

0 

t3 

o 

O 

0 

O  P 

0  p 

0 

<M 

a 

1  t»  ft  G 
I  P 

ft  ft  G 
I  P  O 
G  G 

£  EPS 

O  *H  p 
H  ft  cd 

*  >  I  £ 

O  ft  P 
rH  G 
ft  ft 
ft  G 

§G  p 
G 

t>>  G 

§ 

P  0  0 
hi  tiii) 
t»  G  ft 
O  P  0 
G  ft  P 


GPP 
O  G  G 
G  0  G 
G  CT< 

0  G  G 

rH 

4h  0  Cm 

3  5  ° 


gft  ^ 


•  G  p 
G  0 
O  O  G 
G  G  S3 


ft  ft  G 

co  ft  p 

■H  G 

ft  G  0 
GOG 
co  o  G 

G  I*  G 
•H  O  ft 
ft  0 
G  G  p 
P  ft 
GOp 
G  P 
rH  r-H  G 
G  G 

§  ^ 
0  0 


I  bp  & 

ft  GO 
G  ft  H 
G  ft  G 
B  G  > 
O  CJ  G 

a  g  ft 


w  a 

C5  O 

•  *H 
P  co  p 
0  CJ  O 
p  o  g 

0  -H 
§  -P  0 
CD  O  ^ 
•H  -H  d 
H  H  -P 


cd  co  P 

'S  ^  ft 

G  O 

g  a 

ft  O  G 

p  a  o 
S’gfS 

ft  ft  -rH 
p  ft  ft 
G  G  ft 

ft  3  G 

•H  p  bD 
G  G  ft 
G  ft  G 
O  ft 
O  G 
G  ft 
•'ft  ft 
ft  P 
<U  ft 
G  ft  0 
ft  G  O 

>  G 

G  G 
0  ft  bO 
O  G 
G  rH  ft 

G  ft  P 

>  Ptft 


G  O 
i — I  0  ft 
O  -o  ft 


ft  bD  O 
G  O 
G  ft 

e  8  & 

a  h  o 

p 

G  G  O 
ft  G  G 
P  ft  ft 
O  G 
ft  ft 

°  G  G 
G  ft  G 

O  P 
ft  G 
P  ft  ft 
G  G  ft 
•H  p  > 
P  to  O 
G  G  P 


ft 

o 

I 

— 

1 

•rH 

cm 

C  G  IX! 

i 

•H 

a 

i 

G 

0 

P 

O 

cd  P 

0 

P 

0 

rH  bD 

o 

O 

G 

P 

cd 

G 

P 

P 

p 

CO 

G  G 

p 

•H 

O  p 

cd 

0 

co 

cd 

Ph 

cd 

CO 

•H 

cd  co 

P 

P 

rH 

P  P 

J>> 

G 

CO 

Td 

cd 

P 

0 

Ph 

0  G 

G 

G 

0 

CO 

G 

0 

bD  co 

0 

n 

•s  O 

G 

•H 

P  bD 

O 

P 

G 

y 

G 

CO  *H 

P 

•H 

G  <M 

P 

O 

0  *H 

•H 

P 

1 — 1 

P 

rH 

a 

P  ft  G  ft 
O  P  P 
ft  ft  G  ft 
ft  P 
ft  G  G 
G  ft  ft  S 
G  p  G  P 


C  P  O  P  G 

<U  G  G  ft 

a  G  bp  O  p 

ft  0  ft  O 

P  G  ft  ft 
G  O  ft  ft  O 


O  ft  G 

ft  ft  >>  ft 
p  G  ft  O 
Q  0  ft 

a  g 

•'ft  G 
G  ft  O  O 


ft  0  G 
I  I 

G  ft  G 
P  G  ft 


ft  p  a  o 
X  M  3  d  to 

G  ft  ft  ft  P  ■ 
ft  ‘H  t 

ft  >0  G  G  ft 
O  P  i — t  G  ft 
O  ft  ft  ft  I 

G  p  G  G 

G  O  P  G  p 
ft  G  G  p  FH 
p  ft  ft  G  G 

G  G  ft  p  ft  . 
G  ft  G  G 

p  G  ft 
G~wJ  O  C  G 


G  G 

B?.S£ft 

•HP  G 

Pt  G  G 

a  >  •  g 

G  ft  G  0 
ft  P  O  O 
G  G  G 
ft  ft  G  ft 
G  ft  ft 

G  P  G 
G  ft  G 
G  „  ft 

d  o  a  a 

I  C  O  0  bp 

i  g  a  p  g 

'ft  ®  ft 

l  CT1  G  G  ft 

•  g  ft  a  g 

)  p  0  O  o 
:  ft  p  a  g 


C  >  G 
•H  O  ft  bD 
ft  0  G 
G  G  p  ft  ■ 
P  ftp, 
G  O  P  G 
G  P  ft 


ft  G 

0  >o  < 

a  g 


ft  ft  ft  p 


•  G 

G  G  P 
G  ft  G 


G  G  P 
t>  o  G 


u  P  p 
G  G  o 
ft  G  ft 

ft  a 

G  P  C 
to  O 
ft  0  -H 
O  -i-j  p 
P  ft  G 
p  G  G 
G  C 
O  G  G 

O  ft  FH 

g  a 

p  bD  O 
G  O 
G  FH 

1  §  fr 

)  a  ft  o 
p 

G  C  O 
ft  G  G 
p  ft  ft 
O  G 
ft  ft 
O  p  P 
G  G 
G  ft  G 
O  P 
•H  4) 

P  ft  ft 
D  G  G  -H 
r  -h  ft  > 

I  p  G  O 
i  G  G  p 

)  b»  bD  ft 


•H  ft 
G  ft  O 
ft  ft 
l  ft  G  G 
1  i — I 

P  G  ft 
o  O  ft  -H 
1  ft  ft  CJ 

G 

I  G  p  -H 
'GOP 
O  ft  ft 
i  *H 

CO  Td  0 
H  CO  CJ  x! 

0  cd  -p 
ft  CJ  bD 

■S-SJf 

$  § 
•_a _ 

G  r 

G  G  ft 
■H  t>  G 
ft  G  t> 

•H  G  -H 
CJ  ft  P 


>  P  G 
P 

ft  ft  G 
P  O 
G  G 

t 

O  ft  ft 
ft  ft  G 

.  t  a  > 


ft 

C  G  CO  I 

G  P  G 

G  G  P 
to  G  ft 
P  ft 

G  ft  G  ft 
ft  bD  O  ft 
G  ft  ft  - 

§rH  to  > 
ft  ft  I 
P  G  P 
G  ft  G  G 


O  O  ft  p 

•h  ft  X  bn 

G  G  G 
S  ft 

ft  °  ft  t>, 
FH  O  Sn 
co  O 
ft  G  co  ft 
O  O  G  CJ 
ft  ft  G 
G  P  p  ft 
CJ  G  G  to 
C  0  to  ft 
G  a*  p 
0  G  G  G 

rH  tO 

ft  G  ft 

•5  5  °  5 


£  ft  3  G 


to^  3>. 

ft  G  co 


>  G  to 
G  G  M  CO 
§  ft  G  ft 
0  P  ft 
ft  G  P 
S  ft  H 
G  to  O 
ft  ^  P 
P  >>  P  • 
P  c  to 
ft  G  O  to 
O  G  CJ  G 
^  ft  G 

t>3  a  p  g 

P  ft  G  > 
ft  ft  G  ft 

ft  g  a  p 

ft  P  O  CJ 

ft  cp  a  g 

G  ft 

P  G  ft 
Pt  •  ft  G 
G  co  ft 
od  P  ft 
G  C  ft  O 
GOP 

g  a  P 

ft  O  G  G 
p  a  O  O 

g^p 

•H  »H  *H  C 

^  (O  'h  OJ 

0  U  *H  B 
red  ^3  CJ  O 
*H  -p  bD  S 

CO  CO  *H 
C  *H  CO  0 


0  Xi 

•N  -p  4h 
rH  -p  O 
0  -p 

to  H  d  c 


Td  -p  m  a3 
O  cT  *H 
CO  rH  *H  p 
!  cd  *h  Jh  cd 
1  >  Ph  X  > 


UNCLASSIFIED 


UNCLASSIFIED 


